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Cs3TleCly is shown to possess a rhombohedral unit of 
structure with a9=9.58A, a=83° 48’, containing two 
molecules. The crystals exhibit extensive and quite perfect 
60° twinning on the basal plane, so that most specimens 
simulate complete sixfold symmetry, but when this is 
taken into consideration, a structure based upon 
D;/—R 3 accounts well for observed intensity data. The 


structure may be described conveniently as an aggregate 
of Cs* and T1,Cls~* groups packed together in such manner 
that all large atoms fit into an approximately closest- 
packed arrangement. The Tl2Cl9~* group is considered to 
be a covalent complex. The derived structure is compared 
with the known structures of Cs;As2,Cly and KMgFs. 


INTRODUCTION 


N the course of a search for compounds of 

trivalent thallium which would lend them- 
selves to structural investigations, we have ex- 
amined cesium enneachlordithalliate, CssTl2Clo, 
first prepared by Pratt.1 While this compound 
can scarcely be considered remarkable for tie 
simplicity of its structure, it has proved more 
amenable to treatment than other thallic salts, 
which seem in general to crystallize with highly 
complex structures. We have determined from 
x-ray data the atomic arrangement existing 
within crystals of Cs;TleClo, leading to some 
interesting results with regard to the coordination 
and bonding properties of TI'™ in addition to 
accounting for the stability of crystals of this 
formula. 


DETERMINATION OF THE STRUCTURE 


In accordance with the results of Pratt, we 
were able to prepare clear stable crystals of 
Cs;Tl2Cly with two different habits: Very thin 


'Pratt, Am. J. Sci. 149, 397 (1895). See also Groth, 
Chem. Kryst. Engelmann, Leipzig 1, 437 (1908). 


uniaxial hexagonal plates and hexagonal prisms, 
depending upon the relative concentrations of 
CsCl and TICI; in the solutions from which the 
compound crystallized. The compound is highly 
absorbing, but as the plates were less than 0.1 
mm in thickness, it was possible to get x-ray data 
from various types of transmission photographs. 
We have for the most part worked with the 
plates, although such data as were obtained for 
the prisms are in agreement with the derived 
structure. 

Most of our Laue photographs, prepared with 
the incident beam exactly perpendicular to the 
basal plane, show no deviation from the sym- 
metry of Dey. In two cases, however, we found 
exceedingly thin crystals which gave excellent 
Laue photographs showing only the symmetry of 
Dsa. With either type of crystal the smallest 
hexagonal unit which accounts for the numerous 
data of our Laue and oscillation photographs is 
one with do= 12.8A, co= 18.30A. Our value of the 
axial ratio, co/ao= 1.43, agrees within the limit of 
error with the crystallographic value, Vv 3-0.8257 
= 1.430. 

That we have to deal with extensive and quite 
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perfect twinning on the basal plane is rendered 
highly probable by the following considerations: 
We find that the sixfold symmetry of the one 
type of photograph can be completely reproduced 
both as to the positions and relative intensities of 
all spots by the simple process of superposition of 
two of the threefold patterns after one of them 
has been rotated through 60° about c. We simply 
average the intensities of {HJ-L} and {HI-L} 
of the threefold pattern to get the corresponding 
intensities of either form in the sixfold pattern. 

All of the data of the threefold photographs, 
moreover, are accounted for by the rhombohedral 
unit with one-third the volume of the hexagonal 
unit. This unit, which has ap=9.58+0.04A, 
a= 83° 48415’, contains 2 Cs3TlsCly and leads 
to a calculated density of 4.32 g/cc, is to be 
accepted as the true one. 

The regular vanishings? observed are {hhl} 
with / odd. The space-groups accounting for 
these are C3,°—R 3 c and D3a®—R 3 

The crystal radius* for coordination number 
six of Cl- is 1.81A. A hexagonal unit of 6 layers 
made up of 72 spheres of this size in closest- 
packing would have ap= 12.54A, co=17.72A. Our 
hexagonal cell contains 54 Cl- and 18 Cs*, the 
latter with radius 1.69A only slightly smaller than 
that of Cl-. As we find ap=12.8A, co=18.3A, it 
appears highly probable that cesium and chlorine 
ions are to be found in a nearly closest-packed 
arrangement in Cs3Tl2Cly. Remembering that 
within this hexagonal cell we must have the 
rhombohedral lattice points at 334 and 423, and 
adding the quite reasonable assumption that two 
cesium ions are not to be in contact with each 
other, we are led to a unique distribution 
(excepting small distortions) of the 18 Cst+ and 
54 Cl- within the hexagonal cell. We may 
conveniently symbolize this ideal arrangement, 
which is made up of identical layers, as follows: 
Let A represent a closest-packed layer containing 
3 Cs and 9 Cl with Cs at 02, 34, 20, B a similar 
layer with Cs at 03, 3%, 30. The sequence of 
layers along c is then ABABAB.:- -, correspond- 
ing to hexagonal closest-packing (see Fig. 1). 

If thallium and cesium atoms are not to be in 


2 Astbury and Yardley, Tabulated Data for the Examina- 
tion of the 230 Space-Groups by Homogeneous X-Rays, Phil, 
Trans. 224, 221 (1924). 

3 Linus Pauling, J. Am. Chem. Soc. 49, 765 (1927). 
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Fic. 1. Packing drawing of two layers parallel to (00-1) 
showing the undistorted closest-packed structure. Cesium 
may be distinguished from chlorine atoms by the type of 
shading. 


contact, and the resulting structure is to possess 
rhombohedral symmetry, the 12 TI are to be 
placed in approximate positions given by 
Z2+1/12 and +5/12 relative to each rhombo- 
hedral lattice point. Each Tl is then surrounded 
by 6 Cl. This undistorted or ideal structure can 
be developed from either C3,° or D3a°, in both 
cases with highly specialized parameter values. A 
careful study of the structure reveals no reason 
for eliminating the centers of symmetry and 
twofold axes, and we accordingly assume D3,° to 
be the correct space-group. 

Using the rhombohedral notation of Wyckoff* 
proper to D3a°, we have the following assignment 
of positions: 4 Tl in 4c, uuu, etc. 6 Cs in 6¢, vio, etc. 
6 Cl in 6e, wwo, etc. 12 Cl in 12f, xyz, etc. We ex- 
pect u$ 1/12, y 23, 220. 

The basal plane reflections deserve special 
consideration since they involve but two param- 
eters. Only the even orders appear of course, 
for which 


S (nnn) =4fr1 cos 6rnu+ 6(fostfe,) 
+12fc1 cos 2rn(x+y+2). 
Intensities were calculated from the formula 
1+ cos? 26 
2 sin 20 


The f-values of Pauling and Sherman’ are used 


4R. W. G. Wyckoff, The Analytical Expression of the 
Results of the Theory of Space-Groups, 2nd Ed., p. 157. 
Carnegie Institution of Washington, 1930. 

5 Pauling and Sherman, Zeits. f. Krist. 81, 1 (1932). 
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CRYSTAL STRUCTURE OF Cs;3T12Cls 


throughout this paper. No temperature factor is 
included. The extremely thin plate-like but very 
strongly absorbing crystals used for these 
photographs were so small as always to be 
completely bathed in radiation; the effective 
cross section taken from the beam was pro- 
portional therefore to sin 6, or to n. 

A comparison of calculated with observed 
intensities for the first four orders from (111) 
with u=0.098, x+y+zs=}3 (the value for un- 
distorted closest-packing) is given in Table I. 


TABLE I. Comparison of calculated with observed intensities 
for (nnn) reflections. u=0.098, x+-y+z=}. 


Calculated 
intensity* 


Observed 
intensity 


16 
24 
8 
1 


201 
(nnn) reflections except for the omission of the 
factor n. The intensity comparisons of Table II 
provide strong evidence that we have arrived at 
essentially the correct structural type, although 
small deviations from the ideal parameter values 
are probably indicated. A change of +0.01—0.02 
from the value 3 of the cesium parameter, which 
affects the calculated intensities quite appreci- 
ably, is certainly dimensionally admissible. A 
change in either direction decreases somewhat 
the calculated intensities of all the reflections 
with appreciable structure factors of Table II 
excepting the last three, which are increased. We 
have not been able to decide in which direction 
the parameter value should be moved, so that we 
are giving interatomic distances in Table III for 


TABLE III. Interatomic distances in Cs3T]2Cls. 


* Temperature factor not included. 


The agreement is quite satisfactory, but is 
seriously impaired if the value of u be changed by 
more than about 0.002. A small deviation of the 
chlorine parameter (x-+y+z=3Z, where Z is the 
parameter along the hexagonal axis) from the 
ideal value } is to be expected, although exerting 
a nearly inappreciable effect upon the calculated 
intensities. 

Table II contains a comparison of calculated 
with observed intensities for a group of reflections 
which do not depend upon the thallium parameter 
using the closest-packed values for cesium and 
chlorine. The technique here employed was that 
of transmission through a thin plate which was 
completely bathed in radiation, and the formula 
used in calculation was similar in form to that for 


_ TABLE II. Comparison of calculated with observed intensi- 
for reflections independent of u. x=3, Y=, 
s= w= 


Intensity 
Obs. Calc.* 


Intensity 


* Temperature factor not included. 


Neigh- 
bors 


3 Cl, 
3 Cli 
2 Ch 
2Ch 
2 Clr 
4 Cl, 


Neigh- 
Separation 


3.39A 
3.77 
3.75 
3.80 
3.78 
3.68 
3.68 


Separation 


2.66A 
2.54 
3.70 
3.83 
3.76 
3.62 
3.74 


The intensity data of Tables I and II were 
obtained from twinned crystals, but include only 
planes which possess structure factors inde- 
pendent of the type of twinning postulated. 


DISCUSSION OF THE STRUCTURE 


A diagram illustrating the closest-packing of 
the cesium and chlorine atoms of two layers is 
shown in Fig. 1. Disregarding the thallium 
atoms, which are introduced in positions of 
rhombohedral symmetry, the fundamental trans- 
lation along c of the closest-packed arrangement 
is just the thickness of two layers. As placed in 
the derived structure each Tl is surrounded by 
6 Cl, but pairs of TICl, octahedra share faces to 
give Tl,Cly groups. In the two molecule rhombo- 
hedral unit, we find a Tl.Cl, group of symmetry 
D; at the origin and an equivalent group at $33. 
These groups fit together along the threefold 
axes to form closest-packed columns, which are 
bound together by cesium ions in such manner as 


— 

; 2 16 
| 4 29 
6 18 
8 4 
| 
TI 276 
Tl Cl, 2Cly 
Cs Ci, 2 Clr 

Cs Cl, 2 Chi 
Cs Cla 2 
| Cs 2 Clu 
Cs Cli 1 Che 
| 

(hkl) Obs. Calc.*| (hkl) 

10T 24 24 633 9.7 

202 0 0.04} 13. 

303 16 2.5 

404 0 0.09| 415 10 

505 <1 15 527 | 4.8 

606 6 ©6118 213 5.6 

2iT 48 67 325 1.2 

422 36 54 314 23 


to complete the closest-packing in three di- 
mensions. 

The radius ratio of the univalent radii,® 
r71+3/%co:-= 0.635, indicates that thallic thallium 
should show a coordination number of six (or 
perhaps more) toward chloride ion if the ionic 
model be a good approximation. The mere 
existence in this structure of the complex Tl2Cl, 
groups, in which TICl, octahedra share faces, is 
enough to show that the ionic model is unsatis- 
factory. It seems very improbable that a solution 
containing a mixture of truly ionic Tl**, Cl- and 
Cst+ (or ionic complexes) would give rise to 
crystals in which the thallium atoms are placed so 
near together in pairs, although crystals of this 
composition have actually a rather wide range of 
stability. We conclude therefore that the TleCly 
ion is essentially a covalent complex, and we find 
that an assignment of parameter values for 
chlorine differing but little from the ideal set 
together with the experimentally determined 
value for the thallium parameter give rise to 
interatomic distances in agreement with this 
view. With x=0.336, y=0.176, z= —0.005 and 
w=0.153, we find TI—Cl=2.66A, Tl—Cly 
=2.54A and the CI—Cl separations given in 
Table III. These TI—Cl separations may be 
compared with the sum of the ionic radii,’ 2.76A, 


Linus Pauling, J. Am. Chem. Soc. 49, 765 (1927). 
Ibid. 51, 1010 (1929). 
7 Linus Pauling, J. Am. Chem. Soc. 49, 765 (1927). 
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and the covalent bond distance, 2.55A, obtained 
by adding the “octahedral” radius* of Tl'™ to 
the covalent radius of chlorine (for which the 
“tetrahedral” and ‘‘normal” radii are equal). 
Because of the mutual repulsions of the (some- 
what ionic) thallium atoms within the complex 
and to the less favorable bond angle, 79° 20’ as 
against 96° 20’, the TI—Cl bonds involving the 
shared chlorine atoms are somewhat weaker than 
the others. The value of w is assumed so as to 
make the shared face of the complex have some- 
what shortened edges, in analogy to what is 
observed in truly ionic crystals. A change in this 
value of w together with a possible small change 
from the value of the thallium parameter as 
determined would modify slightly the bond angle 
and interatomic distance for Tl—Cl. 

The structure of Cs3TlsCly is related to the 
structures of Cs3AseCly® and KMgFs; by virtue of 
the fact that all of them are based upon closest- 
packing of the large cations and anions, which 
occur in the ratio one to three. Characteristic 
molecular aggregates are, however, very different 
within the three structures, for we find, respec- 
tively, TleCl>~*, AsCl; and MgF¢~* occurring as 
rather well-defined groups. 

We are indebted to Professor Linus Pauling for 
invaluable suggestions during the course of this 


8 Pauling and Huggins, Zeits. f. Krist. 87, 205 (1934). 
9 Hoard and Goldstein, J. Chem. Phys. 3, 117 (1935). 
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The Crystal Structure of Diaspore* 


F. J. Ew1nc, Gates Chemical Laboratory, California Institute of Technology 
(Received January 21, 1935) 


The parameters in diaspore are redetermined, and the structure is discussed from the 
standpoint of the coordination theory. The existence and location of hydrogen bonds in the 
structure are established, and representative formulas are proposed. 


INTRODUCTION 


HE mineral diaspore, AIHOs, is described as 
forming orthorhombic crystals, with the 
axial ratios a:b :c=0.9372 : 1 : 0.6039, ac- 
cording to Groth. Its composition, symmetry and 
axial ratios suggest its isomorphism with the 
goethite modification of FeHOe, and with 
manganite, MnHOs. A second modification of 
FeHOs, lepidocrocite, also exists and forms 
orthorhombic crystals. The existence of a second 
modification of AlHOz corresponding to lepido- 
crocite is indicated by the x-ray studies made by 
Bohm! on certain bauxites. 

Work begun by the writer in 1929 had disclosed 
the structure of diaspore and suggested a 
probable structure for lepidocrocite. The com- 
pletion of this work was unavoidably delayed 
until recently, and in the meantime the structure 
of diaspore was derived by Deflandre? and 
independently by Takané,? and Goldsztaub‘ 
found that the goethite structure was essentially 
that of diaspore. 

In the present paper the diaspore-goethite 
structure is discussed from the standpoint of the 
coordination theory, with particular reference to 
the distribution of the hydrogen atoms and their 
relation to the existence of a second modification.® 
The results of a more accurate determination of 
the parameters in diaspore are given, and the 
interatomic distances thus determined are shown 
to establish the existence and location of the 
hydrogen bonds as suggested by the coordination 
theory. 


* No. 456. 

1J. Bohm, Zeits. f. anorg. Chemie 149, 203 (1925). 

* M. Deflandre, Bull. Soc. Franc. Mineral. 55, 140 (1932). 

*Takané, Proc. Imp. Acad. Tokyo 9, 113 (1933). 

*S. Goldsztaub, Comptes rendus 195, 964 (1932). 

* The structure of the second modification as typified by 
lepidocrocite will be described in a future publication. 


TuHeE UNIT oF STRUCTURE, SPACE GROUP SyYM- 
METRY AND THE GENERAL POSITIONS 
OF THE ATOMS 


The unit cell of diaspore was determined from 
oscillation and Laue photographs to correspond 
to the axial lengths 


a=4.40A, b=9.39A, c= 2.84A. 


These values give the axial ratios 0.469 :1: 
0.302, in good agreement with the crystallographic 
ratios 0.4686 : 1 : 0.3020 obtained by doubling 
the b-axis. The observed density of diaspore, 3.4, 
indicates that a unit of these dimensions contains 
4 AlIHOs. 

The space group of diaspore was determined 
from Laue data to be V;,!*—Pbnm, as found by 
Deflandre and Takané. 

The observed even orders of reflection from 
(001) were found to decline normally in intensity, 
indicating that all the atoms are situated in 
planes perpendicular to c at a distance c/2 apart. 
To place the oxygens anywhere except on the 
reflection plane m would give O—O distances of 
1.42A or less, which is impossible. All the atoms 
therefore lie on the m reflection planes, the 
general positions, referred to a center of sym- 
metry as origin, being as follows: 


4 Alatu,v,—4;3—u,3+2, —4;3+u, $—v, 457, 0, }. 
40, at “1, 01, —}; etc. 
4 Oy at U2, v2, —}; etc. 


THE DETERMINATION OF PARAMETERS 


Oscillation photographs were taken of Mo Ka 
radiation by reflection from (010), ten even 
orders of reflection being observed. Using the 
usual methods it was found that the observed 
intensities of reflections at large angles serve to 
fix v, the aluminum parameter, at +0.144+0.003 
irrespective of the oxygen positions. By using 
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values of the aluminum parameter within these 
limits it was then found possible on the basis of 
observed intensities to fix the oxygen parameters 
within the ranges —0.216<v1<—0.190 and 
— 0.062 <ve< —0.042. The parameters were then 
adjusted within these limits to give the best 
quantitative relation between observed and 
calculated intensities. The values thus obtained 
are 


v= +0.146+0.001, —0.200+0.002, 
and ve= —0.052+0.002, 


in which the probable limits of error as given 
correspond to the minimum displacement causing 
disagreement between observed and calculated 
intensities. The extent of the agreement between 
observed and calculated intensities of reflection 
from (010) using these values of the parameters is 
shown in Table I. 


TABLE I. Observed and calculated intensities of reflection on 
oscillation photograph from (010) with [001] as axis. 


Intensity 


Intensity 
s.* Calc. 


Obs.* Calc. 


5.4 2.0 


i) 


tnt 
on 
SSSLS 


7.3 
4.7 
1.2 
2.1 
0.1 
0.4 
0.0. 
0.0 
0.1 


COM 


(0.20.0) 


* Estimated visually, and corrected in the case of side spectra for 
varying time of reflection. 
b Based on Bragg and West's f values. 


The relative shortness of the a-axis limits the 
amount of diffraction data obtainable by re- 
flection from (100), so that an exact determi- 
nation of the parameters along a by rigorous 
methods is difficult. By assumptions which will 
be discussed in the next section it is possible to 
predict values of these parameters, the predicted 
values being 


u= —0.036, ui= +0.27, 


These values give excellent quantitative 
agreement between the calculated and observed 
intensities of reflection, as is shown in Table II. 
A shift of 0.01 in any of these parameters some- 
what impairs the quantitative nature of the 
agreement. 


—0.21. 


EWING 


TABLE II. Observed and calculated intensities of reflection on 
oscillation photograph from (100) with [001 ] vertical. 


Intensity Intensity 
Calc.» Obs.* Calc.» 


® Estimated visually, and corrected in the case of side spectra for 
varying time of reflection. 
b Based on Bragg and West's f values. 


The parameter values determined by Deflandre 
are, restated in our convention, w= —0.042, 
+0.25, —0.25, v= +0.144, —0.20, 
ve= —0.05. His values for v, 7; and v2 are in close 
accord with ours, but his values for u, u; and 12 
do not account quantitatively for reflections from 
(100) as well as do ours. This is illustrated, for 
instance, by the reflection (701), which is 
observed 0.6 as strong as (501); the calculated 
ratio is 0.5 using our parameters values and 0.2 
using Deflandre’s values. Takané’s parameter 
values differ more widely from ours, and lead in 
some instances to calculated intensities that are 
in qualitative disagreement with our data. 


DISCUSSION OF THE STRUCTURE 


The diaspore-goethite structure can be regarded 
as built of aluminum or iron centered octahedra 
arranged in the manner shown in Fig. 1. The 
octahedra share edges in such a way as to form 
one-dimensionally infinite aggregates which may 
be called ‘double rutile strings” running parallel 
to the c axis. These strings are joined by sharing 
corners, as shown in the figure. Each anion is 
common to three octahedra. 

This type of structure is readily predictable on 
the basis of the coordination theory® if a distri- 
bution of the hydrogens is assumed giving all the 
anions the form (OH;)~3, such anions being 
required by the electrostatic valence rule to be 
shared between three octahedra.? The only 
alternative distribution of the hydrogens gives 
the anions O= and OH-, and the electrostatic 


6 Linus Pauling, J. Am. Chem. Soc. 51, 1010 (1929). 

7 With the character of the polyhedra and their method of 
sharing defined, the number of structures fitting the unit 
cell is greatly limited, and the structure shown in Fig. 1 1s 
the only one that is immediately apparent under the condi- 
tions given. 


(200) 0.2 0.13 ie 0.8 0.81 
(400) 5.0 3.4 | (301) 1.7 2.2 
(600) 0.0 0.00 | (501) 0.5 0.57 
(500) 0.0 0.08 | (701) 0.3 0.28 
(020) 4.0 
(040) 10 S 
(060) 6.0 (061) 
(080) 0.25 (081) 
(0.10.0) 3.0 (0.10.1) ‘ 
(0.12.0) 0.7 (0.12.1) 
(0.14.0) 1.1 (0.14.1) 
(0.16.0) 0.0 (0.16.1) ‘ 
(0.18.0) 0.25 (0.18.1) 
0.7 (0.20.1) t 
2 
a 
0) 


CRYSTAL STRUCTURE OF DIASPORE 


A 
q 
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Fic. 1. The ideal structure of diaspore, showing the 
“double rutile strings’ formed by the aluminum-centered 
oxygen octahedra. The relationship between the octahedra 
and the close packed arrangement of the oxygens is also 
shown. The aluminums are represented by the small 
spheres. 


valence rule requires that O= and OH~ anions be 
shared between four and two octahedra, re- 
spectively. This latter type of sharing is actually 
that found in the structure of the second modifi- 
cation of FeHOs, lepidocrocite.* 

The coordination theory not only gives the 
ideal structure successfully, but also indicates the 
degree of distortion to be expected. Edges shared 
between octahedra, according to the rule for 
distortion, should be reduced to 2.50A, with 
other edges being compensatorily lengthened. 
The lengthened edge is given in the case of 
diaspore by the c identity distance, 2.84A. By 
fitting the distorted octahedra to the unit cell 
values for the entire set of parameters can be 
obtained, and v, v; and ve thus obtained are 
within 0.02 of the values derived rigorously from 
intensity measurements. 

Starting with the rigorously determined values 
of v; and ve the requirement that shared edges be 
2.50A long leads directly to the values of +0.27 
and —0.21 for u; and we, respectively. As previ- 
ously shown arguments based on quantitative 
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agreement with intensities serve to confirm these 
values, but equally important is the fact that 
they introduce into the structure the 2.50A 
shared edge which appears in many accurately 
known structures such as rutile, anatase, brookite, 
and corundum, and which has been shown to 
correspond to a minimum of crystal energy in 
the case of rutile and anatase.*® 

The value of -—0.036 for the aluminum 
parameter u is calculated by making the reason- 
able assumption that structurally similar Al—O 
distances inside the octahedron are equal, i.e. 
that the aluminum is equidistant from the three 
oxygens in the center of the double rutile string, 
and likewise equidistant from the three oxygens 
at the corner of the string. These two conditions 
lead to slightly different values for u, the average 
of which is given by —0.036. 

The formal existence of (OH,)~? anions is most 
readily explained by postulating hydrogen bonds 
between pairs of oxygens. Hydrogen bonds, 
forming collinear OHO groups, have been found 
in certain crystals containing oxygen and acidic 
hydrogen, namely H;BO3;,9 NaHCO;,!° and 
KeHPO,," and in and the crystalline 
hydrates 2H2O and LisSO,- The 
hydrogen bond has also been shown to exist in 
gaseous (HCOOH). by means of electron diffrac- 
tion methods.'® 

In diaspore the location of the hydrogen bonds 
can be expected to be between oxygens in 
different octahedra because of the effect of 
Al—H repulsion; their exact location is finally 
determined by a comparison of the interatomic 
distances concerned. 

Interatomic distances calculated from our 
values of the parameters are shown in Table III. 
With the exception of the shared edges (which 
are 2.50A long), the O—O distances within each 
octahedron lie in the range from 2.73 to 2.84A, or 
somewhat in excess of 2.70A, the sum of the ionic 
radii. Between different octahedra the O—O 
distances are 3.10A or greater with the important 


8 Linus Pauling, Zeits. f. Krist. 67, 377 (1928). 

9 W. H. Zachariasen, Zeits. f. Krist. 88, 150 (1934). 
10 W. H. Zachariasen, J. Chem. Phys. 1, 634 (1933). 
J. West, Zeits. f. Krist. 74, 306 (1930). 

1929), H. Barnes, Proc. Roy. Soc. (London) 125, 670 
13 W. H. Zachariasen, Zeits. f. Krist. 89, 442 (1934). 
4G, E. Ziegler, Zeits. f. Krist. 89, 456 (1934). 

% L. Pauling and L. O. Brockway, Proc. Nat. Acad. Sci. 

20, 336 (1934). 
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TABLE III. Interatomic distances in diaspore. 


Neighbors in Neighbors in 


Atom same octahedron other octahedra 
O; 2 O; at 2.84A 2 On at 4.11A 
4 O; at 2.78 2 On at 3.92 
2 On at 2.78 1 On at 2.71 
2 Oy at 2.73 
1 O;, at 2.50 
On 2 O,, at 2.84A 20, at 4.11A 
2 O, at 2.78 20, at 3.92 
2 O,; at 2.73 2 O, at 3.10 
2 Oy at 2.50 10, at 2.71 
1 O,; at 2.50 
Al 1 O, at 2.01A 2 Al at 3.10A 
2 On at 1.98 
1 O, at 1.86 
2 O, at 1.82 


exception of the 2.71A distance. This distance of 
2.71A is not only characterized by the fact that 
it is much the shortest of similar distances, but 
also by the fact that it is inexplicably so unless 
it is assumed to correspond with a hydrogen 
bond. 

The location of the 2.71A distances in the 
diaspore structure is shown by the small tubes in 
Fig. 2. This figure perhaps makes it clear that by 
slightly twisting the double strings about their 
points of contact the distances corresponding to 
the tubes can be lengthened and brought into 
equality with other O—O distances between 
octahedra. The only apparent explanation for the 
structure having a twist which differentiates 
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these distances is that the tubes correspond to 
hydrogen bonds."® 

The allocation of one hydrogen to each tube 
gives the correct stoichiometric composition to 
the crystal, and hence accounts for all of the 
hydrogens. 

The value of 2.71A found for the OHO distance 
in diaspore is in agreement with the values 2.71 
and 2.67A found in H3BO3; and (HCOOH)., 
respectively, and lies somewhat above the values 
2.54 and 2.55A found in KeHPO, and NaHCO;, 
and below the values 2.77 and 2.96A found in 
H.O and LisSO,-H.O. This indicates that the 
strength of the hydrogen bond in diaspore is 
considerably greater than that in ice or in the 
hydrate LieSO,-H.O, and is approximately the 
same as in and (HCOOH). 

The Al—O distances in diaspore correspond 
closely to the Al—O distances of 1.99 and 1.85A 
found in corundum, since in diaspore the 
Al—Oy,, distances average 1.99A and the Al—O, 
distances average 1.84A. The slight deviations of 
0.01 to 0.02A from these average Al—O distances 
shown in Table III are within the error of the 
parameter determinations, but on the other hand 
they are in the proper direction to be related to 
a distortion produced by the hydrogen bond. In 


16 This differentiation of distances is not seriously af- 
fected by variation of the parameter values over ranges 
considerably greater than their probable error. Deflandre’s 
values of the parameters might have been used to deduce 
the hydrogen end, obtaining an OHO distance of 2.60A 
in comparison with a next smallest distance of 2.80A. 


Fic. 2. A pictorial representation of the diaspore structure, with the hydrogen bonds shown as tubes. 
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both corundum and diaspore the shift of the 
Al—O distances from 1.90A, the sum of the radii 
concerned, is evidence of the effect of Al—Al 
repulsion. 

The formula for diaspore we have written 
AlHOs, the HOsz group so indicated being taken 
as representing hydrogen-bonded oxygens. The 
formula HAlO2 would be equally correct as 
indicating that both H and Al cations formed 
coordinated structures with the oxygen cations, 
and the chemical inference in this latter formula 
that diaspore is meta-aluminic acid is borne out 
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by the existence of meta-aluminate salts such as 
the spinels. Either formula would be more 
correctly written in the polymerized form 
(AIHO2), or (HAIO2), to admit the known 
coordination number of 6 for Al. The usual 
formulas for diaspore, AlxO3;-H2O, AlxOo(OH)s, 
or AIO(OH) are wrong in indicating a distinct 
existence of water or OH™ ions. 

I am gratefully indebted to Professor Linus 
Pauling for his original suggestion as to the 
probable structure of diaspore and for his 
interest and advice during the work. 
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On the Normal Vibration Frequencies of NH;, PH; and AsH; 


J. B. Howarn, Chemical Laboratory, Harvard University 
(Received February 9, 1935) 


It is shown that the assumption of valence forces for 
NHs3, PH; and AsH,; leads to the conclusion that the high 
frequency, degenerate normal vibration, although the- 
oretically active, does not appear as a fundamental in the 
infrared spectra of these compounds. This conclusion 
follows from the fact that the remaining three funda- 
mentals, all of which have been observed, can be accounted 
for satisfactorily with proper values of the two constants 
of the potential function for valence forces. The values 
calculated for this absent frequency with the same force- 
constants are 3450 cm, 2340 cm™ and 2125 cm™ for 
NH3;, PH; and AsHs, respectively. The frequency for NH; 
may be fixed more precisely at 3470 cm™ if we interpret 


HE justification for using the concept of 

valence forces in treating NH; is based 
upon its reasonably successful application to 
H.O and CHy. J. H. Van Vleck and P. C. 
Cross! have shown that the fundamental fre- 
quencies of HO can be computed with fair 
accuracy using a valence-type potential function. 
Jenny E. Rosenthal? has found that the fre- 
quencies of CH, are consistent with such a 
potential function to within six percent. It is 
therefore justifiable to make the assumption 
that a similar kind of forces exists in NH; and 


1983 H. Van Vleck and P. C. Cross, J. Chem. Phys. 1, 357 
"ies E. Rosenthal, Phys. Rev. 46, 730 (1934). 


the 4417 cm™ band as a combination band involving this 
frequency and the 950 cm™ frequency. Other evidence is 
given in support of a value for the frequency in the neigh- 
borhood of 3450 cm~. There is, however, the difficulty of 
explaining why the predicted frequencies are not ob- 
served. The 4417 cm™ band is probably to be ascribed to 
the fourth fundamental in case this fundamental does 
actually appear. If this is the case, then there is an extreme 
departure from valence forces in NH; and probably likewise 
in PH; and AsH;. Frequencies are computed for NH2D, 
ND:H and ND; to provide a means of testing the value 
predicted upon the assumption of valence forces for the 
fourth fundamental in NHs3. 


that the valence-type potential function is 
probably a good first approximation to the 
general quadratic expression. Any large deviation 
from this approximation would, on this assump- 
tion, then be considered a consequence of the 
anharmonic terms introduced by the passing of 
the nitrogen atom through the plane of the 
hydrogen atoms. Consequently, to afford a com- 
parison with the results obtained for NH3, we 
have also considered PH; and AsHs3, in which 
the bonds are doubtless similar to those of NH; 
and in which this so-called tunnel effect is 
thought to be of less consequence.* 


3See for example D. M. Yost and T. F. Anderson, J. 
Chem. Phys. 2, 624 (1934). 
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CALCULATION OF FREQUENCIES 


Under the assumption of valence forces the 
potential function of the symmetrical pyramidal 
molecule YX; is 


V=1/2K(AR?+AR2?+ AR;’) 
+1/2HR,?(Aais? + Aaes?+ 


Here K and H are force constants, AR; is the 
departure of the distance Y—X; from its equi- 
librium value Rp and Aa;; is the departure of 
the angle X -Y—X; from its equilibrium value ao. 
This molecule has four distinct frequencies: »; 
and v3, in the notation of Dennison,‘ denote the 
high and low frequencies, respectively, of vibra- 
tions with change of electric moment parallel to 
the figure axis. Likewise v2 and v4 correspond to 
the high and low frequency vibrations, re- 
spectively, with change of electric moment per- 
pendicular to this axis, the latter vibrations 
being doubly degenerate. All four frequencies 
are theoretically active in both the Raman and 
infrared spectra. 

Formulas for calculating the fundamental fre- 
quencies of vibration have been derived by 
Lechner® and others for the above potential 
function. In Table I are given the values of K 
and H which make the calculated values of »; 
and v4 agree with the experimental values we 
have assigned to these frequencies. The values 
computed for v2 and v3 with these force constants 
are likewise given; in the case of NH; this calcu- 
lated value of v3 is the unstarred frequency. 
It is seen that the difference between the calcu- 
lated value of v3; and the observed frequency in 
its neighborhood, which we have assigned to 73, 
is rather large for NH; but is smaller for PH; 
and AsH3. 

The large discrepancy in the case of NH; can 
be explained as a consequence of anharmonic 
terms in the potential function and does not 
otherwise, invalidate the assumption of valence 
forces. These terms principally affect the v3 (950) 
frequency and necessitate a large correction in 
obtaining the observed frequency, which corre- 
sponds to a vibration with finite amplitude, from 
the frequency for infinitesimal amplitude of 
vibration, called the mechanical frequency. Mr. 


4D. M. Dennison, Rev. Mod. Phys. 3, 280 (1931). 
5F. Lechner, Wien. Berichte 141, 633 (1932). 


TABLE I. Bond angle ao,° observed infrared frequencies’ 
and calculated frequencies in cm, and calculated force con- 
stants in dynes/cm for NH;, PH; and AsH3. 


vs (||) 


(1) 
v1 
ve (L) 


1121 
2327 


3428 


(1005) 
(2122) 
2125 


Kx10% 
HX 10-5 


6.35 
0.59 


* Corrected for anharmonic term. 


M. F. Manning has kindly consented to calculate 
the mechanical frequency directly from his 
symmetric double-minima potential function® 
for NHs3. He finds it to be 1050 cm-, in contrast 
to the experimental frequency of 950 cm™. 
The value 1110 cm, calculated for »; upon the 
assumption of valence forces, does not, of course, 
allow for the anharmonic correction. It will 
resemble the mechanical frequency since the 
anharmonic corrections for the other frequencies 
are no doubt much smaller than that for 73. 
Applying Manning’s correction to the computed 
value of v3 we obtain a final calculated frequency 
of 1010 cm™, which is to be compared with the 
experimental frequency, 950 cm-'. Thus the 
agreement between the observed frequency 7; 
and the frequency calculated upon the basis of 
valence forces is as good for NH; as for PH; and 
AsH;, in which there is probably no large 
anharmonic correction necessary. Furthermore 
the agreement is to within about six percent. 
This gives one a certain amount of confidence in 
the valence-force approximation. 

No frequencies have been observed in the 
infrared, however, in the neighborhoods of the 
values calculated for ve, excepting those assigned 
to 

6 The angle for NH; is known from spectroscopic data. 
The angles for PH; and AsH; were chosen on the basis of 
theoretical considerations and evidence from crystal struc- 
ture given by Pauling, J. Am. Chem. Soc. 53, 1367 (1931). 
No appreciable change is produced in the force constants 
by a change of a few degrees in the angles. 

7 Source of frequencies: NH;—Schaefer and Matossi, Das 
Ultrarote Spektrum, p. 250. PH;—Fung and Barker, Phys. 
Rev. 45, 238 (1934). AsH;—Robertson and Fox, Proc. Roy. 


Soc. A120, 149 (1928). 
8M. F. Manning, J. Chem. Phys. 3, 136 (1935). 
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DISCUSSION OF v2 
NH; 


With the value of v2 calculated on the assump- 
tion of valence forces it is possible to explain 
the 4417 cm infrared band as the combination 
ve (3450) +3 (950) =4400 and, indeed, it 
is difficult to account for the band with any 
other combination of frequencies. The structure 
of the band is very irregular and does not admit 
a simple analysis.* Because of this complexity it 
is probably not a parallel-type band. Such bands 
have, in general, a relatively simple structure 
and are free of the distortions resulting from 
rotation-vibration interaction. Hence we must 
account for what is probably either a_per- 
pendicular-type structure or an overlapping of 
several bands. There appear to be only four 
alternative explanations. (a) The assignment 
(3335)+ 3 (950)=4285 cm suggested by 
Badger and Mecke!® leads to a parallel-type 
band and would not alone explain the structure. 
Now it is possible that the second overtone of 
vs (1631), viz., 3v,=4893 (minus the anhar- 
monic terms), is stronger than the first overtone 
and accordingly might badly distort the parallel- 
type structure of this combination band. This 
alternative is improbable, however, because of the 
wide separation of the two frequencies involved. 
(b) There is the possibility that the second over- 
tone of v4 gives rise by itself to the 4417 cm™ 
band although this is doubtful since the band 
is almost as intense as the fundamental at 3335 
cm-!. (c) The combination v2 (3450)+ 73 (950) 
= 4400 cm involving the valence-force value of 
ve does, on the other hand, account for the fre- 
quency and intensity of the band, and leads to 
a perpendicular-type structure, the complexity 
of which would be increased by the doublet 
character introduced by the 950 cm~! component. 
(d) The possibility that the 4417 cm~ frequency 
is itself v2 is considered below. 

In the Raman spectrum of gaseous NH; 
Lewis and Houston" have observed a line at a 
shift of 3219 cm— which they suggest may be re. 
The difference between this frequency and the 


°E. F. Barker, private communication. 

1°R. M. Badger and R. Mecke, Zeits. f. physik. Chemie 
BS, 333 (1929). 
(1935) M. Lewis and W. V. Houston, Phys. Rev. 44, 903 
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valence-force value calculated for ve might con- 
ceivably be ascribed to the presence of cross 
terms of moderate size in the potential function. 
There are two bits of evidence, however, which 
conflict with the allotment of this line to vp. 
Firstly, Bhagavantam™ finds that the corre- 
sponding line with 3208 cm™ shift in the Raman 
spectrum of the liquid is very well polarized, a 
result which is incompatible with the assignment 
of this line to a degenerate vibration. Secondly, 
the combination v2 (3219)+ 3 (950) =4169 
gives too low a frequency for the 4417 cm~ band 
which is difficult to account for by another 
assignment. 

Further support of the valence-force assign- 
ment for v2 is given by the frequencies of the CH; 
group in the methyl halides. The CH;X mole- 
cule has two more frequencies than the NH;- 
type molecule as a result of the additional 
atom X. If we allow the bond between C and X 
to vanish, however, these two frequencies will go 
to zero and, simultaneously, the vibrations of 
CH;X will pass over into the vibrations of a 
free CH; molecule. In the spectra™ of the methyl 
halides the four frequencies characteristic of the 
CH; group are readily identified. The strength 
of the C—X bond in each halide can be estimated 
from the observed frequency of vibration of the 
CH; group against X by using the diatomic 
model (CH;) — X. For X equal to F, Cl, Br and I, 
the force constant for stretching of the C—X 
bond assumes the values 5.3X10°, 3.3105, 
2.8X10° and 2.2X10*° dynes/cm, respectively. 
Plotting these values against the observed fre- 
quencies of the CH;X molecules and extrapo- 
lating to zero force-constant, we obtain the 
following approximate frequencies of vibration 
of a free CH; molecule: »3=900 cm-, »4= 1420 
cm-!, 2985 cm and ve=3120 These 
values are similar to the frequencies calculated 
for NH; upon the assumption of valence forces. 
The extrapolation does not, of course, take into 
account the consequences of the tunnel effect. 

The most obvious difficulty with the valence- 
force value of 3450 cm™ for v2 is that such a 
frequency does not appear in the infrared whereas 


2S. Bhagavantam, Ind. J. Phys. 5, 49 (1930). For a 
discussion of the significance of the polarization of Raman 
lines see G. Placzek, Leipziger Vortrage, p. 71 (1931). 

18 The frequencies of the methyl halides are given by A. 
Adel and E. F. Barker, J. Chem. Phys. 2, 627 (1934), 
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ve is a theoretically active frequency. Very little 
is known, however, regarding the change of 
electric moment during vibration in such mole- 
cules as NH; and it cannot be said how serious 
a difficulty this is. The failure to appear in the 
Raman spectrum of the gas may be due merely 
to insufficient intensity ; v4 (1631) is likewise not 
observed. It is commonly true that theoretically 
active, degenerate frequencies appear less intense 
in Raman spectra than frequencies of completely 
symmetric vibrations and in the case of com- 
pounds containing hydrogen are frequently not 
observed. 

Another possible difficulty has been pointed 
out to the author by Professor Dennison. From 
the complete expressions derived by Salant and 
Rosenthal" for the symmetrical tetratomic mole- 
cule Y°X;™ we obtain for two such molecules 
(denoted by single and double primes) the rela- 
tion = where 
b= 6u(c?/a?) +1 and w= M/(3m+ M). Here c and 
a are the altitude of the pyramid and the side 
of the base, respectively. This relation must be 
satisfied rigorously as long as the potential 
functions of the two molecules remain purely 
quadratic. Assuming that the 4417 cm~ band 
of NH; and the 3287 cm~ band of ND;" arise 
from the combination v2+v3, we obtain 3470 
cm and 2520 as the values of v2 for 
and NDs, respectively. These frequencies, when 
substituted in the above relation together with 
the observed values of v4, lead to a ratio of 
0.94 instead of 0.87 which is the value of the 
ratio of the (b/u) term for NH; to that for NDs. 
This difficulty is serious only if the centers 
ascribed to the observed bands for v2+v3 and 
v4 are correct and if our method of obtaining 
values for v2 is valid. It does not require a large 
change in the values used to yield a ratio 0.87. 
Thus the values ve=2565 cm-! and »4=1180 
cm™ calculated below for ND; assuming valence 


4 J. E. Rosenthal, Phys. Rev. 45, 426 (1934); 47, 235 
(1935). These formulas are given in terms of the force 
constants of the general potential function. In making the 
calculation, the general constants were evaluated in terms 
of K and H by means of the expressions for the former given 
by E. O. Salant and J. E. Rosenthal, Phys. Rev. 42, 812 
(1932) on p. 819. To obtain agreement with the valence- 
force formulas’ for YX; it was found necessary to change 
the negative sign before Ke in the expression for D to a 
positive sign. 

1 Observed ND; frequencies taken from Silverman and 
Sanderson, Phys. Rev. 44, 1032 (1933). 
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forces give the correct ratio (since they: are ob- 
tained from the NH; frequencies by using the 
formulas from which the above comparison 
relation is derived). 

The 4417 cm~ and 3287 cm frequencies, on 
the other hand, give (as was pointed out by 
Professor Dennison) a ratio of 0.89, in closer 
agreement. A check of this sort does not, of 
course, imply that the frequencies used are the 
correct ones since only the ratios of the fre- 
quencies are tested. Nevertheless, this bit of 
evidence, together with the difficulty of explain- 
ing the 4417 cm~! band as a combination other 
than v2 (3470)+ 3 (950)=4420 cm, indicates 
that if the assumption of valence forces does not 
give the correct value of ve, this frequency of 
4417 cm™ is probably the fourth fundamental. 
If such is the case then there is a very large 
departure from valence forces in passing from 
CH, or H2O to NHs3. 


Pp 'H; and AsH; 


No evidence is available regarding the valence- 
force interpretation of ve in PH; and AsHs. It 
seems probable from structural considerations, 
however, that the values of v2 in Table I, cal- 
culated using valence forces, will be equally good 
for these molecules and for NH3. 


CALCULATION OF THE ISOTOPE EFFECT 


If the valence-force frequency, ve= 3450 cm“, 
for NH; is approximately correct the assumption 
of valence forces should be sufficient to calculate 
the fundamental frequencies of the isotopic 
molecules NH2D, ND2H and ND3. In making 
this calculation we have used the complete 
formulas derived by J. E. Rosenthal" for the vi- 
brational frequencies of the symmetric tetratomic 
molecule YX2X*. Values of the force-constants, 
K=6.42X10° and H=0.59 X 10°, were chosen so 
as to give the calculated frequencies v2= 3470 
and »v4= 1631 for NH3. This value of 
ve satisfies the relation v2+950 cm—!= 4420 cm™. 
Table II gives the results of the calculations. 
The letter A after a frequency of NH2D or 
ND:H indicates that it corresponds to a vibra- 
tion which is antisymmetric with respect to 
reflection through a plane passing through the 
figure axis and the hydrogen isotope which 
differs from the other two. The other frequencies 
of NH2D and ND-H are symmetric with respect 
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TABLE II. Observed" and calculated fundamental frequencies 
in cm™ of NH; and isotopic molecules. 


3335 ? 
3355 (3470) 


1631 
(1631) 


obs. 950 
calc. 

1110 
1025 1385A 1600 3400 2510 3470A 
925 1230 1460A 2455 2565A 3435 
770 1158 2421 ? 


840 1180 2405 2565 


NH; 

NH.D_ obs. 
calc. 
obs. 


calc. 


to this plane. The computed values of v3 will lie 
above the observed frequencies, as is found to be 
the case for NHs3, and will more closely resemble 
the frequencies of vibrations with infinitesimal 
amplitude. 

The agreement between the frequencies calcu- 
lated for ND; and the three observed frequencies 
is satisfactory. As we have noted above the band 
reported at 3287 cm= may be used to determine 
ve and gives the value »x2= 2520 cm™, in rough 
agreement with the calculated value. 


OVERTONE-COMBINATION BANDS 


The proximity of »; and v2 in all three mole- 
cules would explain the fair success of efforts to 
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interpret the spectra using only the three funda- 
mentals which appear. An attempt to obtain 
further evidence regarding the valence-force 
value of ve from overtone-combination bands is 
rather inconclusive as long as the structures of 
these bands or the actual appearance of the 
bands to be expected from certain combinations 
of fundamentals are not better known. Thus, 
for example, we have two possible explanations 
of the 5059 cm= band of NH3. The assignment 
given by Badger and Mecke!® is », (3335) 
(1631) =4966 cm“. The alternative assign- 
ment is ve (3470)+ 4 (1631)=5101 cm. The 
former combination leads to a perpendicular- 
type band, the latter to a band with both 
perpendicular and parallel components. The 
structure of the band has been analyzed as per- 
pendicular by Stinchcomb and Barker’ but it 
is possible that the presence of a parallel com- 
ponent would be masked by the more complex 
perpendicular component. 

The writer is indebted to Dr. E. Bright 
Wilson, Jr., and Professor D. M. Dennison for 
their helpful criticisms and suggestions in con- 
nection with this work. 


16 Stinchcomb and Barker, Phys. Rev. 33, 305 (1929). 
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The Magnetic Susceptibility of Lithium Hydride 
The Diamagnetism of Ionic Crystals 


SIMON FREED AND Harry G. THODE, George Herbert Jones Laboratory, University of Chicago 
(Received January 7, 1935) 


The magnetic susceptibility of crystalline lithium hydride was measured at room temperature 
(300°A) and at 78°A and was found to be —4.60X10~® per mole with practically no tem- 
perature coefficient. One theoretical prediction deviates from the experimental value by 90 
percent and another by 130 percent. This discrepancy is interpreted as due mainly to a para- 
magnetism independent of the temperature which is present in all crystals and which has not 


been included in the calculations. 


HE diamagnetic susceptibility of a mole of 
atoms or atomic ions may be expressed by 


x= 


where x is the susceptibility per mole, N is 
Avogadro’s number, 7? is the mean square 
distance of the 7’th electron from the nucleus, 
the summation extending over all electrons and 
e, m, c have their usual meanings. Several theo- 
retical methods have been devised for evaluating 
72; one is based on Hartree’s method of self- 
consistent fields,! another is due to Van Vleck 
and Pauling? independently, and more recently, 
Brindley? has employed the approximate eigen- 
functions of Slater. Angus? has also computed r? 
with these eigenfunctions but in greater detail. 
The diamagnetic susceptibilities derived with 
these various methods are compared by Angus 
with the susceptibilities of the alkali halides in 
solution as measured by Ikenmeyer.? The sus- 
ceptibilities of Pauling are from 10 percent to 
60 percent and of Stoner from 50 percent to 80 
percent higher than the experimental values. 
Most of the values of Brindley are also high 
while those of Angus fluctuate above and below 
the experimental values by a few percent. 

All the theoretical methods deal with electronic 
systems subject to central forces such as isolated 
ions in field-free space and since the only data 
available refer to ions in crystals or in solutions, 
there is ever present the question of the validity 
of such comparisons. The chief reassurance has 
come from the (supposed) additivity of the ionic 
susceptibilities, suggesting magnetic independ- 
ence between them, and further there seems as 

1Stoner, cited by W. R. Angus, Proc. Roy. Soc. A136, 


569 (1932). 
2 Cited by W. R. Angus, see reference 1. 


yet little reason to suspect considerable changes 
in the mean square radius when a gaseous ion 
takes its place in a crystal lattice. 

In the hope of clarifying some of the un- 
certainties we measured the magnetic suscepti- 
bility of lithium hydride principally because it is 
the simplest of all ionic lattices.* And when its 
eigenfunction is determined a direct comparison 
with experiment will be available. It remains to 
be seen whether the eigenfunction of the crystal 
will be separable into the eigenfunctions of the 
ions in the lattice with sufficient accuracy. This 
separability is implicit in all previous investiga- 
tions, a matter we shall discuss later. In terms 
of separate ions, it may be observed that the 
hydride ion is probably more sensitive than any 
other to changes in electrical environment. Hence 
a measurable difference in the susceptibility 
would probably accompany a change from the 
gaseous state to the crystalline and the ion might 
serve to indicate whether all ions undergo such 
changes. 

MATERIALS 


We have had the good fortune to obtain 
crystals of pure lithium hydride. They were 
furnished by the Maywood Chemical Works, 
Maywood, N. J., through the courtesy of Dr. 
H. Osborg. In the light of chemical and magnetic 
analyses, any ordinary paramagnetic impurity 
was present in less than 0.002 percent, probably 
much less, and diamagnetic impurities such as 
lithium hydroxide, sodium hydride made up less 


3 Hylleraas (Zeits. f. Physik 63, 771 (1930)) has confirmed 
from the heat of sublimation that the lattice is ionic. The 
crystals are similar to sodium chloride, colorless, hard, 
melt at a rather high temperature (680°C), all in keeping 
with an ionic lattice. It was found to be cubic, NaCl type, 
(Bijvoet and Karssen, Proc. Amst. Acad. 31, 496 (1922)). 
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MAGNETIC SUSCEPTIBILITY OF LITHIUM HYDRIDE 


than 0.5 percent of the material and could be 
ignored in magnetic considerations. Chemical 
analyses for the ferromagnetic impurities, iron 
and nickel, were made as these more than others 
were likely to contaminate the material on 
account of the method of manufacture. The 
presence of total iron, both free and combined, 
was found to be less than 7 moles iron in 10 
million moles lithium hydride and probably less 
than 1 in 10 million. No trace of nickel was 
observed although 1 mole nickel in 5 million 
lithium hydride would have been detected. It is 
well known that ferromagnetic metals commonly 
lose their ferromagnetism when alloyed so that 
the best sign of ferrornagnetic contamination is 
the magnetic behavior of the substance, the 
necessary criterion being a dependence of the 
susceptibility on the field strength. We are unable 
to observe such a dependence in our materials. 

After the lithium hydride had been measured 
magnetically, it was taken out of the tube for 
chemical analysis. Different samples gave 99.4 
percent, 97.5 percent, 99.5 percent an average of 
98.8 percent lithium hydride. The impurities are 
probably lithium hydroxide and other alkali 
hydrides. The hydroxide is formed from the 
hydride and moisture. For this reason the utmost 
precautions were taken never to expose the 
hydride to moisture. All the necessary manipula- 
tions (such as grinding the crystals to a powder, 
tamping it into a glass tube for the magnetic 
measurements) were carried out in a ‘“‘dry-box”’ 
through which dried gases were flowing. In 
selecting the crystals, any having a blue colora- 
tion were avoided. This color is a very sensitive 
indication of uncombined lithium metal. It is 
probable that even if some of the faintly blue 
crystals had been taken no difference would 
have entered into the magnetic results within 
our limits of error. 

About six months after the first magnetic 
measurements were taken, we received some 
lithium hydride from the Maywood Chemical 
Works which had been made in a more recent lot. 
Dr. Osborg was good enough to have crystals 
selected for us which had little likelihood of 
having touched the metallic walls of the appa- 
ratus in which it was made. From these we 
further selected colorless crystals for the mag- 
netic measurements. 


MAGNETIC MEASUREMENTS 


The well-known method of Gouy was em- 
ployed for the magnetic measurements and at 
low temperatures the method described by 
Freed. The magnet was of the Weiss type 
kindly loaned us by Professor W. D. Harkins. 
It was built with the aid of a grant given to 
Professor Harkins by the Joseph Henry fund of 
the National Academy of Sciences. The first two 
sets of runs listed below were carried out with 
magnetic fields up to 23,000 gauss and in the 
remaining runs with fields up to 18,000 gauss. 

The magnetic susceptibility was also deter- 
mined at 78°A to enhance the effect of any para- 
magnetic impurity since paramagnetic suscepti- 
bility usually varies almost inversely as the 
absolute temperature. That is, the paramagnetic 
impurity would have its room temperature sus- 
ceptibility almost quadrupled at 78°A. 

As has been noted, we were unable to detect 
any variation in the susceptibility with field 
strength such as one would expect from the 
presence of any ferromagnetic impurity. 

On May 21, 1934, the magnetic susceptibility 
per mole of lithium hydride was determined at 
room temperature (300°A) as —4.4210~°. 

On May 22, 1934, another sample of the same 
crystals was taken and independent magnetic 
measurements at room temperature (300°A) 
gave —4.68 x 10-* per mole. The deviations from 
the average of these two sets of runs is 3 percent, 
a little high, perhaps because of the irregularities 
in the packing of the powder in the tubes. 

On December 10, 1934, lithium hydride which 
had been prepared afresh by the Maywood 
Chemical Works was measured under conditions 
which were very different from those which had 
prevailed in the other measurements. The glass 
tube containing the powder was changed, the 
pole gap of the magnet and hence range of mag- 
netic fields were different. The susceptibility per 
mole was found in a series of runs as — 4.65 X 10-® 
in satisfactory agreement with the values ob- 
tained six months previously. These crystals had 
been measured in a Dewar vessel at room tem- 
perature. The temperature was then lowered to 
78°A (the temperature of liquid nitrogen) and 
the susceptibility at this temperature was — 4.60 


4S. Freed, J. Am. Chem. Soc. 52, 2702 (1930>. 
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Fic. 1. Change in apparent weight in magnetic field. Water for calibration: upper curve, weight at left. 
Lithium hydride: lower curve, weight at right. Clear circles at room temperature 300°A. Halved circles at 


78°A. 


<10-*, an increase of 1.1 percent. This difference 
is real as the experimental method is especially 
suitable for measuring differences resulting from 
changes in temperature. At room temperature, 
the percentage of the susceptibility which may 
be ascribed to a paramagnetic impurity is about 
(78/222) K1.1=0.35 percent. The deviations 
between measurements within a set of runs 
occurred in the third decimal place. 

The figure gives the curves obtained from the 
measurements both of water, which served for 
calibration, and of lithium hydride at 300°A 
and at 78°A. 

We conclude, then, that the molal suscepti- 
bility of lithium hydride in the crystalline state 
is —4.6010-* within about 2 percent. 


DISCUSSION OF RESULTS 


The molal susceptibility of lithium hydride 
predicted by Pauling is —8.65X10-* and by 
both Brindley and Angus — 10.54 10-6 and the 
experimental value is —4.60 x 10~*, truly another 
order of magnitude. We are faced with several 
possibilities of which the more important are: 
(1) The mean square radii of the electrons in the 


ions (hydride ion, especially) have undergone a 
shrinkage of 50 percent or more; (2) a factor 
which does not operate in isolated ions introduces 
in crystals a paramagnetism which is practically 
independent of the temperature. On this basis 
the measured diamagnetic susceptibility cannot 
be related exclusively to the term of the mean 
square radii. This factor, we are inclined to 
believe, accounts for most of the discrepancy.® 


5 That there may also be a possible shrinkage in the 7; of 
the electrons is suggested by the fact that the radius of the 
hydride ion as defined for x-ray data is less in lithium 
hydride than in other alkali hydrides. Zintl and Harder 
(Zeits. f. physik. Chemie B14, 265 (1931)) give the follow- 
ing lattice distances for the alkali hydrides. Upon subtract- 
ing the radii of the positive ions as given by Zachariasen 
(Zeits. f. Krist. 80, 137 (1931)) from one-half the lattice 
distance, the radii of the hydride ions are obtained. These 
are listed in the last column. 


Lattice 
distance 


It need not follow that there is any difference in the r;? of 
the electrons in the various hydride ions. Even if the effect 
on the 7;? were as large as is here suggested, it would be 
small compared with the discrepancy between the com- 
puted and experimental susceptibilities and we are inclined 
- the view that the ignored possibility (2) is a dominant 
actor. 
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MAGNETIC SUSCEPTIBILITY OF LITHIUM HYDRIDE 


Van Vleck® has shown that in atomic systems 
having spherical symmetry, the instantaneous 
angular momentum of the 4S») state is zero 
as well as its average angular momentum. 
In polyatomic systems, among which we may 
include crystals, the instantaneous angular mo- 
mentum is not zero, even when the stationary 
state is 1X» that is, even when the average 
angular momentum is zero. The perturbation 
theory foresees a paramagnetic contribution from 
the square of these instantaneous fluctuations in 
angular momentum which of necessity is positive. 
If the frequency differences between the various 
states is high as compared with kT this para- 
magnetism would be independent of the tempera- 
ture. The magnitude of this ‘‘fluctuation’’ para- 
magnetism which we shall call the “Van Vleck 
paramagnetism” is in many instances great 
enough to account for the discrepancy between 
the experimental susceptibility and the various 
computed values. As examples we may cite the 
magnetic measurements’? on KMnQ,, 
KeMoO, and UO2SO, which were selected to 
show this effect in MnO,-, Cr2O;-, MoOg and 
UO,**. In these instances the paramagnetism 
was more than great enough to overbalance the 
total diamagnetism of the crystal and yet it 
was small compared with the usual paramag- 
netism which varies with the temperature. 

There is no doubt that the alkali halides may 
undergo quantum transitions* to states associ- 
ated with angular momentum. The same is 
doubtless true in lithium hydride’ also. If we 
symbolize the NaCl lattice as 

NatCl-Nat 
Cl-NatCl- 
Na*Cl-Nat*- 


upon absorbing ultraviolet light, the crystal is 
activated into a state thus, 


NatCl-Nat 
Cl-Na°Cl° 
NatCl-Nat 


to which the blue color has been ascribed.® Other 
transitions are possible also. 


*J. H. Van Vleck, Electric and Magnetic Susceptibilities, 
Oxford University Press, 1932. 
— and C. Kasper, J. Am. Chem. Soc. 52, 4671 
* R. Hilsch and R. W. Pohl, Zeits. f. Physik 68, 721 (1931). 
11932) Bach and R. F. Bonhoeffer, Naturwiss. 20, 940 
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The final state we have been discussing is 
paramagnetic but it is not the number in the 
final state we are concerned with. It is rather 
the fluctuation from the basic state to the 
activated state which accounts for the para- 
magnetism independent of the temperature. The 
conclusion is then that some of the diamagnetism 
originating in the mean square radii is masked 
by the Van Vleck paramagnetism in all ionic 
crystals as well as in others; and that it is 
erroneous to employ experimental diamagnetic 
susceptibilities alone to furnish average radii of 
ions or to serve as criteria for the accuracy of 
such approximate eigenfunctions as have been 
proposed. In some instances, it is true such 
errors may not be great, a matter which requires 
separate investigation. Theories which furnish 
correct mean square radii would yield greater 
diamagnetic susceptibilities than are to be found 
experimentally and hence the values of Pauling 
or of Stoner should not be rejected because of 
“poor agreement” with the experimental sus- 
ceptibilities. 

From this point of view, we would expect the 
additivity of ionic susceptibilities not to hold 
except as a rough approximation especially 
among dissimilar compounds. One would expect 
some difference between the diamagnetism of 
crystals and of their solutions and possibly also, 
but to a much lesser degree, some departure 
from linearity in the susceptibility as a function 
of the concentration in regions of high con- 
centration. 

Very opportunely, Hoare” has remeasured 
with considerable care the diamagnetic suscepti- 
bilities of salts with rare gas ions and we shall 
quote some of his conclusions : 


It will be seen then that no conclusive evidence as to the 
validity of the additivity principle can be derived from 
this table, especially from Ikenmeyer’s results; the results 
of the present investigation indicate, however, that it may 
be true to a first approximation. 

Hocart measured the susceptibilities of the chlorides in 
the crystalline form and in aqueous solution, finding the 
results in the former state to be less than the latter. It 
would appear, therefore, from the results for the bromides 
and iodides, that the result found by Hocart for two salts 
might very well be in agreement with a general principle. 

This conclusion is strengthened by Kido’s results. 


10F, E. Hoare, Proc. Roy. Soc. A147, 88 (1934). 
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A Plant for the Production of Heavy Water 


W. G. Brown anp A. F. DaGGett, Department of Chemistry, Columbia University 
(Received February 14, 1935) 


A description is given of a plant for the production of deuterium oxide by the electrolytic 


method which operates with a higher efficiency than has been previously reported for this 
method. The fractionation factor, uncorrected for evaporation, averages 8.6. Steel cells, which 


also serve as the cathodes, are used and the entire construction is explosion-proof. 


EVERAL types of laboratory equipment for 
the production of deuterium oxide by the 
electrolytic method have been described.' To 
these published accounts we wish to add the 
following brief description of a plant recently 
constructed in this laboratory which may be of 
interest because we are able to report a higher 
efficiency than is usually realized in the electro- 
lytic separation of hydrogen and deuterium. 
This plant consists of 37 water-cooled steel 
cells of 3 1 capacity each, operated in series at 100 
amperes, and arranged in three stages, and one 
smaller cell of 800 cc capacity which constitutes 
the fourth and final stage. The gases generated in 
each stage are collected by iron header pipes and 


Net 


LLL 


Fic. 1. Cross-sectional diagram of cell. 


1G. N. Lewis and R. T. MacDonald, J. Chem. Phys. 1, 
341 (1933); H. S. Taylor and H. Eyring, Proc. Am. Phil. 
Soc. 72, 255 (1933); H. S. Taylor, H. yg ee A. A. 
Frost, J. Chem. Phys. 1, 823 (1933); W. D. Harkins and 


C. Doede, J. Am. Chem. Soc. 55, 4330 (1933); L. C. Ander- 
son, os Halford and J. R. Bates, J. Chem. Phys. 2, 342 
(1934. 
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are conducted successively through spray traps 
and explosion traps into air-cooled burners in 
which recombination takes place. The entire 
construction is sufficiently rugged to withstand 
explosions without damage. 

It was intended that this plant should operate 
on water already enriched to the extent of 0.5 
percent deuterium oxide and with this as the 
starting material four stages suffice to increase 
the concentration to 99.9 percent, the daily 
production then being approximately 90 grams. 
However, the inadequate supply of 0.5 percent 
water frequently necessitates operation on water 
of considerably lower concentration and the 
average production falls far short of the figure 
mentioned. 


DESIGN OF THE CELLS* 


The cells are constructed of Shelby seamless 
steel tubing, 4’’ in diameter and 18”’ in length. 
A bottom plate and a central tie-rod which 
secures the lid are attached by arc welding. The 
lid carries a gas outlet tube and a tube through 
which the cell is filled and drained (these tubes 
are not shown in the accompanying sketch). 
The cell itself, insulated from the cooling water 
by paint, serves as the cathode. The anode con- 
sists of a perforated cylinder of sheet nickel, 
34” in diameter and 5” high, attached by nickel 
rivets to the flattened end of a 3” nickel rod 
which passes through an inverted rubber stopper 
in the lid of the cell. The nickel cylinder is 
insulated from the cell wall by small pieces of 
hard rubber. The small cell which constitutes 
the fourth stage is constructed on the same 
design but the material is 23” iron pipe. 

The current density at the cathode varies with 
the level of the liquid in the cell. At the anode 


*The design of the cells was worked out by Dr. D. 
Rittenberg and Dr. S. H. Manian. 
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HEAVY WATER PLANT 


the current density is approximately 0.28 am- 
pere/cm? when this electrode is completely 
immersed, and it increases to about twice this 
value at the end of a run when only half 
immersed. 

Before use the cells were cleaned by elec- 
trolysis of a hot solution containing 30 grams 
each of borax and trisodium phosphate per 
liter. This treatment together with the care 
taken to ensure the purity of the solutions used 
subsequently in the cells appears to be effective 
in preventing foaming. 


CONSTRUCTION OF THE BURNERS 


The gases generated in the cells of the first 
stage pass into iron header pipes which are 
connected to the cells by short sections of heavy 
wall rubber tubing. There are three separate 
header assemblies for this stage, each collecting 
the gas from ten cells. From these the gases 
pass through spray traps packed with steel wool, 
then through explosion traps filled with 60-mesh 
sand, and finally through jets at which re- 
combination takes place. The first stage burner 
consists of a fan-cooled 3” iron pipe mounted at 
a slight incline. A cap fitting over the upper end 
carries the three converging jets which project 
into the pipe. At the lower end a nipple conducts 
the steam into vertical water-cooled condensers. 
The jets consist of }” iron plugs through which 
holes of 0.090’ diameter have been drilled and 
which are threaded into the 3” iron pipes pro- 
jecting through the burner cap. The explosion 
traps are mounted rigidly on the burner cap and 
are connected to the spray traps by rubber 
tubing thus permitting the cap to be withdrawn 
for lighting. The entire assembly, with the 
exceptions mentioned, is constructed of standard 
plumbing fittings. 

The second stage burner is similar in size and 
construction and the principle of multiple con- 
verging jets is retained by having two jets, each 
fed by three cells, with separate headers and 
traps. The orifice in these jets is 0.050” in 
diameter. 

The third and fourth stage burners are smaller 
in size, each having a single gas jet. The orifices 
are 0.028” and 0.015’, respectively, the sizes 
being somewhat critical with the smaller gas 
flow in these burners. 


ELECTRICAL EQUIPMENT 


The cells are connected in series and are 
operated from a 110-volt line, the current being 
controlled by a rheostat of 0.4 ohm resistance. 
An underload circuit breaker prevents the re- 
sumption of the electrolysis should an inter- 
ruption in the current sufficient to extinguish 
the burners occur. The fourth stage cell takes 
approximately one-half of the current used in the 
other cells (100 amperes), the remaining half 
being passed through a shunt. This stage and 
the third stage are provided with short-circuiting 
switches allowing these stages to be switched in 
or out as desired. The fourth stage is operated 
only when attended, but the remainder of the 
plant is operated overnight regularly without 
attention. 


OPERATION OF THE PLANT 


The thirty cells of the first stage are charged 
with 3 | each of a 3 percent sodium hydroxide 
solution made up of the water to be concentrated. 
The sodium hydroxide is a c.p. Electrolytic 
grade, and the water is distilled once or twice 
before use. When it appears desirable a distilla- 
tion from alkaline permanganate is carried out. 
The charges are measured out and placed in the 
cells with the aid of a portable apparatus built 
for this purpose. 

When the electrolysis has proceeded for ap- 
proximately 72 hours, this charge, totalling 90 
liters, will have been reduced to about one- 
sixth of its original volume. It is then removed, 
a part of it distilled and then transferred to 
the second stage. Similarly the residue from the 
second stage is transferred to the third, and 
so on. A charge of 90 | of water containing 0.5 
percent deuterium oxide will yield 15 | of 
approximately 2.5 percent. In the second stage 
the concentration will have increased to about 
12 percent, in the third to about 60 percent and 
in the fourth to better than 99 percent. 

The disposal of the recombined water depends 
upon its concentration. Normally the water 
used in the first stage contains from 0.3 percent 
to 0.5 percent deuterium oxide and except for 
the first fraction the recombined water contains 
an appreciable amount of deuterium. It is our 
practice to discard recombined water containing 
less than 0.05 percent deuterium, and to re- 


3 
217 

ps 
in 
ire 
nd 
ite 
he 
ise 
ily | 
as. 
nt 
ter 
he 
ire 
th. 
ch 
he 
gh 
eS 
h). 
ter 
yn- 
el, 
el 
od 
er 

is 
of 
tes 
me 
ith 
de | 
D. 


218 W. G. BROWN AND A. F. DAGGETT 


work the remainder. Recombined water from 
the later stages is returned to the previous 
stage. 

FRACTIONATION FACTORS 


Samples of recombined water are collected 
from each stage at the end of each run. These, 
as well as samples from the residues, are analyzed 
interferometrically? for their deuterium content 
providing data from which the fractionation 
factors may be calculated. Data for a number of 
typical runs are given in Table I, together with 


TABLE I. 


Fraction- 
ation 
factor 


First stage 


essssss 


Second stage 


DW OWN 


— 


NOUN 
tn Go Com mein 


1 
2 
3 
+ 
5 
7 
9 
2 
3 


ninth 


Third stage 


the fractionation factors. These are the in- 
stantaneous factors, and in a number of cases 
they have been checked by calculations of the 
overall factors, that is, the factor which must be 
substituted in the integrated Rayleigh distilla- 
tion formula to satisfy this relation between the 
initial and final volumes and concentrations. In 
most cases, however, we have not sufficient data 
to calculate the overall factor since the runs are 


2R. H. Crist, G. M. Murphy and H. C. Urey, J. Chem. 
Phys. 2, 112 (1934). 


rarely straight runs. It happens almost always 
that the electrolysis is stopped one or more 
times during the course of a run to add re- 
combined or other water of appropriate con- 
centration. Thus a first stage run might be 
started with a charge of 0.15 percent water, to 
which would later ‘be added 0.30 percent water, 
and perhaps later 0.50 percent water. 

The average value of the fractionation factors 
listed in Table I is 8.6. This is appreciably higher 
than any previously reported.!:* These values 
have not been corrected for evaporation. Since 
the temperature in the cells averages about 25°C 
such a correction would increase the average 
fractionation factor to approximately 12. Pre- 
sumably this high factor could be realized by 
cooling the gases leaving the cells so as to collect 
water carried over as water vapor, but we have 
not yet attempted to do this. 


PURIFICATION OF HEAVY WATER 


The residue from the fourth stage, after having 
been neutralized and distilled, is not sufficiently 
pure for many experimental purposes. It con- 
tains traces of organic impurities not readily 
removable by distillation or by refluxing with 
oxidizing agents, and to effect the purification 
we make use of an additional stage of electrolysis. 
In this apparatus the deuterium and oxygen 
generated by electrolysis in a small iron cell are 
mixed with an excess of pure oxygen and burned 
at a silica jet inserted in a Pyrex condenser. 
The apparatus is flushed with oxygen before 
electrolysis is begun and with this precaution 
the oxides of nitrogen, customarily present in 
recombined water, are excluded. 

This plant was constructed and is being 
operated under a grant-in-aid of heavy water 
research by the Rockefeller Foundation. In the 
design and construction of it we have enjoyed 
the cooperation of Professor H. C. Urey, Pro- 
fessor J. E. Zanetti, Dr. D. Rittenberg and other 
members of the staff of this laboratory. 


3R. P. Bell and J. H. Wolfenden, Nature 133, 25 (1934); 
B. Topley and H. Eyring, Nature 133, 292 (1934), J. Am. 
Chem. Soc. 55, 5058 (1933), J. Chem. Phys. 2, 217 (1934); 
P. Harteck, Proc. Phys. Soc. 46, 277 (1934). 
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The Thermal Conductivity of Polyatomic Gases 


A. R. UBBELOHDE, Department of Thermodynamics, Clarendon Laboratory, Oxford 
(Received December 12, 1934) 


Reasons are given for supposing that the thermal conductivity of polyatomic gases should not 
be independent of pressure. A calculation is made for carbon dioxide, and the results suggest 
that a portion of the specific heat at 0°C has a relaxation constant considerably less than 10 if 
the discrepancy between experimental and theoretical values of thermal conductivity is due to 
restrictions on the energy exchange bet ween molecules. When such effects are present, determin- 
ations using parallel plates give a more reliable result than those using a vertical wire. 


KNOWLEDGE of the mechanism of heat 
transport in gases is important both in the 
kinetic theory, and in many chemical reactions 
involving activated molecules. The object of the 
present communication is to discuss the transport 
of vibrational heat in gases, and to point out that 
the thermal conductivity of polyatomic gases 
may vary with pressure. This is because of the 
fact that at higher pressures any departure from 
the Maxwell Boltzmann distribution of energies 
is more rapidly readjusted. A calculation is given 
for the heat transport between two parallel 
plates, and when certain simplifying assumptions 
are made, it is shown that the’ thermal con- 
ductivity K at any pressure is related to the 
maximum value Kymax with full participation of 
the vibrational energy, by an equation of the 
form K(1+ F)=Kwmax. The factor F depends on 
the pressure and dimensions of the apparatus, 
and has been evaluated in a typical case for COs. 
The bearing of this calculation on previously 
obtained experimental data for the thermal 
conductivity of polyatomic gases is discussed, and 
it is pointed out that empirical pressure cor- 
rections may be misleading. A correlation of the 
relaxation time for redistribution of vibrational 
energy in COs, as calculated from the dispersion 
of sound at high frequencies, with that calculated 
from thermal conductivity, is shown to be 
reasonable, though a quantitative check is 
hampered by lack of suitable experimental data. 
Two methods are available for determining 
how far vibrational energy of the molecules of a 
gas takes part in heat transport. According to 
views developed by Eucken! if K is the thermal 
conductivity of a gas, C the specific heat at 
constant volume, and 7 the viscosity, the factor f 


‘Eucken, Physik. Zeits. 14, 324 (1913). 


in the equation K=fnC may be calculated. The 
translational energy is assumed to be transported 
according to the Chapman-Enskog formula 
Kerans= 2.59Ctrans, and if both rotational and 
vibrational transport are independent of trans- 
lational, 


K=(Kerans+ Krot+Kviv) 
= (2.5Ctrans+ Crot + Cvin)- 
If y=C,/C, C=R/(y-1), 
and 
Ff (wucken) = (2.5C trans + Crot+ Cviv)/ 
(Cerans t+ Crott+ C vin) 
i.e., f (eucken) =1+1.5Ctrans(¥ — 1)/R=3(97—5). 


If vibration takes no part in heat transport but 
is operative in experiments determining C and y, 


no vib) >= (2.5C trans + Crot) / 
(Crranst+ Crot+ Cyiv) (21/4)(y 1) 


for polyatomic and (19/4)(y—1) for linear 
molecules. 

A test of this calculation requires a precise 
knowledge of thermal conductivity, viscosity 
and specific heats, all at the same temperature 
and pressure. Such a calculation is possible for 
CO:.? For other gases included in Table I less 


TABLE I. 


Gas CHs CeHs CoHe H2S 


1.68 1.57 1.59 1.50 1.72 1.70 
1.69 1.43 1.42 1.54 (1.37) (1. rg 
1.60 1.34 1.37 1.16 1.68 1.63 


(Eucken) 
(no vib) 


2 Kannuluik and Martin, Proc. Roy. Soc. Al42, 362 
(1933). 
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reliable data are available.* Except in the last 
two cases, where there may be serious error in the 
data, the value of f calculated from K/nC lies 
between the two theoretical limits. For (C2H;)2O 
and gases with smaller y the effect is still more 
marked but the data less reliable. Eucken has 
suggested that the value of f is due to the fact 
that vibrational and in some cases even rotational 
energy is not taking full part in heat transport, 
owing to incomplete temperature equilibrium. A 
second method is now available to test this, since 
experiments on the velocity of sound in such 
gases as COz show that at ordinary temperatures 
and pressures dispersion sets in in the neighbor- 
hood of 10° cycles/sec. One explanation is that at 
still higher frequencies the vibrational energies 
are no longer able to follow the adiabatic temper- 
ature changes in a sound wave. 

During the interval of 10-° sec., the average 
distance diffused by a molecule may be calculated 
from the formula #?=2Dt. For COz at 0°C and 
20 cm of mercury, D=0.4 and x=3X10-* cm, 
i.e., some 100 mean free paths. In such a case it is 
clear that in the parallel plate method for 
determining the thermal conductivity of a gas, 
we cannot assume complete equilibrium for the 
vibrational energy in various layers. In order to 
clear up certain difficulties in interpreting 
experiments on the dispersion of sound‘ and on 
thermal conductivity® it is important to correlate 
the two phenomena as far as possible. 

Consider an experiment where the heat trans- 
port in a gas between two parallel plates is 
measured. If Q is the heat transport/unit area, 
and the plates are a distance 29 apart, and are 
maintained at temperatures Tp, and Tp,, the 
experimentally determined thermal conductivity 
is obtained from the equation Q= KAT p/z» where 
ATp=Tp,—Tp,, and K is the mean value of the 
thermal conductivity over small ranges of 
temperature. This value of K may differ from 
that calculated from the kinetic theory of gases 
for one of two reasons: 

(1) The average temperature 7y character- 
izing the distribution of vibrational energies in 
any layer will differ from the temperature T 


3 Landolt Bornstein Tabellen 5th ed. 

4 Richardson, Proc. Roy. Soc. A146, 56 (1933); Richards 
and Reid, J. Chem. Phys. 2, 194 (1934). 

5M. Trautz, Ann. d. Physik 18, 833 (1933). 
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characterizing the distribution of rotational and 
translational energies by an amount L=7T—Ty 
(see below). 

(2) Because of temperature jumps at the 
surface, the difference AT=7,—T>. and still 
more so ATy=Ty,—Ty, may be appreciably less 
than the measured difference AT p= Tp, —Tp,. 

At sufficiently high pressures the problem is 
best considered as a case of diffusion in the 
steady state. If the molecules are distributed 
N,, No/unit volume, in classes characterized by 
vibrational energies F;, Ey (in many cases it is 
sufficient to consider only two levels), then in the 
steady state the rate of increase of N;/unit 
volume by diffusion is 


= aiss (a) 


and must be balanced by the net rate of removal 
within the layer by collisions, i.e., 


aite+(0N1/0t) con =0, 


where (0N,/dt)cou is a function of the compo- 
sition of the layer. For the rotational and trans- 
lational degrees of freedom, (AN,/dt)cou is 
usually assumed to be negligible, because of the 
displacement from a Maxwell Boltzmann distri- 
bution being so small, the use of a temperature 
parameter 7 to characterize this distribution is 
thus allowable. With respect to these degrees of 
freedom we then have for each class 1, 2, etc. 


D(@N,/dz?)=0 and dN,/dz=const=ky,. 


Since the net number of molecules crossing any 
layer in unit time is D(dN,/dz) the heat transfer 


Q==DE,(dN,/dz) =DZE,(dN,/dz) 


if on the average D is independent of class (i.e., 
neglecting the Chapman correction). Integrating 
a, over the distance 29 between the plates, 


If oquilibciem is fully established at the plates, 


so that 


The standard formula is thus obtained. 
We are, however, concerned with a degree of 
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THERMAL CONDUCTIVITY 


freedom for which (0N;/0t)cou is not negligible, 
on account of the finite lag L. The vibrational 
temperature Ty will be defined by Ni/No 
=e~(4:-Bo/k? and the calculation be made for a 
gas where the rate of relaxation is given by 


dTy/dt=B(T—Ty).° (b) 


Furthermore, over the small intervals of temper- 
ature which have to be used in thermal con- 
ductivity experiments 


pCvin= (E,— Ep) (dN,/dTy) (c) 


will be appreciably constant, where Cyj, is the 
vibrational specific heat/unit mass, and p is the 
density. 

In the steady state, the heat transport Q 
through any layer is constant, i.e., 


—Q/NmD=C'(dT/dz)+Cyin(dTy/dz)=b (d) 
where 0 is independent of z, so that 
C'T+ CyinTy =b2+d (e) 


where d is a constant. In (d) and (e), if we use the 
normal coefficient of diffusion, D, C’=2.5Ctrans 
+Crot in view of the Chapman correction. Com- 
bining (a) with (c), and using (b), 


D(d?.N,/dz*) =[pCvin/(E1 — Eo) aist 
Eo) y/8t) con 
= —[eCviv/(Ei—Eo) 8(T—Tv).  (f) 
Substituting for T by means of (e), and noting 
Eo) 
D(@?Ty/dz*) = —(B/C')((b2+d) — (Cvin +C’) Ty). 
The solution to this equation is 


Ty = Ae*?+ (b2+d)/(Crvin +C’), 


where 


(g1) 


a? =B(Cyin +C’)/C'D, 
or since L=T—Ty, from (e) the solution 
simplifies to 


L=Ae**+Be-, (g2) 


where the arbitrary constants are different from 
those in (g;). The third term in (g;) is equivalent 
to a linear temperature gradient for Ty, and the 


® Rutgers, Ann. d. Physik 16, 351 (1933). 
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first two take account of relaxation effects, as 
measured by the constant 8 in the equation 
for 

In order to apply this solution to a concrete 
case, we have to evaluate the arbitrary constants 
by considering the boundary conditions. Since 
the system is symmetrical, when 


z=0, L=L,=A+B 
and when 


L=-L; = Ae**0+ 


whence 


B=-—Ae*™, and L,;=A(1—e***).  (h) 


The constant A can be evaluated by considering 
the change in energy of the molecules by collision 
with the plates. If the lag is appreciable, the 
vibrational temperature of the molecules ap- 
proaching the plate will be appreciably different 
from that of the plate itself, even at ordinary 
pressures. According to the theory of the 
accommodation coefficient’ if Ty be the (vi- 
brational) temperature of the molecules ap- 
proaching the plate, and 7) that of the molecules 
leaving it, 


Tp —Ty =ayin(Tp—Ty), 


where a,j, is the accommodation coefficient for 
vibrational energy. Little is known of the values 
of ayi», but unless the material of the plates 
absorbs the molecules strongly, a,i, is probably 
very small because of the small duration of 
contact with the plate. By neglecting it altogether, 
the maximum value L; can have will be obtained 
and can be calculated. The rate of increase of 
N, ina layer of unit area and thickness dz next to 
the walls will be 


Eo) JD(dTy/dz)o 
+ — Eo) JBLidz 
+(Né/4)[pCviv/(E1— Eo) Javin(T p— Ty) =0. 


In this equation the first term represents the net 
increase of N, by diffusion, which for this layer is 
equal to the net number of molecules crossing its 
boundary, D(dN,/dz)o. The second represents the 
change in J; by collisions in the gas, and the 
third by collisions with the plate. Né/4 is the 


7 Knudsen, Ann. d. Physik 34, 593 (1911), 
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number of molecules approaching the plate/unit 
time. If avi, is negligible, (dT y/dz))>=0. Differ- 
entiating (gi) and writing 


z=0, (dTy/dz).=0, 
(dT /dz))= a(A — B) = aA (1+e**) 
and from (d), making the same substitution, 
C' (dT /dz)o>= —Q/NmD. 


Using these equations, and substituting for A 
from (h), 


(i) 

On inserting limits in (e), 
C'’AT+ CyivATy = Q20/NmD, 

and since 
L2=AT—21i, 

AT(C’+Cyiv) —2L1Cvin=Q20/NmD. 
Eliminating ZL, from (i) and (j), 

AT(C’+Cyiv) =20Q(1+F)/NmD,  (k) 

where 

Since equation (k) may be written, 

Q=(Cviv + C’)(AT/20)-NmD/(1+ F), 
Kmax= NmD(Cyin+C’), 

K=K1ax/(1+F). 


and 


we have 


The only other correction, for the difference 
ATp—aAT is proportional to \/%9 where is the 
mean free path, and is usually negligible at the 
pressures used. 

For a concrete case, inspection of the above 
correction term F shows: 

(1) When a1, indicating a large time of 
relaxation, 


F= (2Cyin/azoC’) (az/2) = 
In this case 
= NmDC’ 


as would be expected when vibration takes no 
part in heat transport. 
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(2) When a is very large, F=0, as would be 
expected when vibration takes its full share in 
heat transport. It is important to know the 
region of pressures over which K is changing with 
pressure. Consider the case of COz2 at 0°C in an 
apparatus for which z9>=0.3 cm 


B=8.1X10° at 30° and 760 mm, 
Cvin= Crr =5.1= C 


By using the data of Richards and Reid‘ for the 
activation energy of vibrational transfer, B= 2.7 
10° at 0° and 760 mm and at a pressure p mm 
Hg, 6=2.7pX 105/760. 


D at 0° and 760mm=0.1 and at p mm Hg, 
D=76/p. 


Since Crr=5 the molecule is taken as linear in 
calculating C’. From these data, 


a? = (2.7 105/76) p? X 11.2/9.5 =4.2p? 


and 


when p=0.1 mm, 1 mm, 10 mm and 100 mm, the 
values of F are 0.152, 0.148, 0.05 and 0.005, 
respectively. The experimental value of F 
calculated from 


1+ F=f (eucken)/K/C 
is approximately 0.05 for COs. 


DISCUSSION 


A comparison of these figures with the 
experimental value 0.05. for CO: which is ob- 
tained in the neighborhood of 200 mm, suggests 
that the value of 8 from acoustical experiments 
may have been overestimated. This is in part 
explained by assuming with Eucken and Miicke® 
that the vibration with @=1830 is not excited 
even at acoustical frequencies of 10° cycles. This 
would lead to a maximum of 0.02 to F at 0°C over 
the whole range of pressures considered above. 

At present a more exact comparison is not 
possible, owing to the fact that the experimental 
methods commonly used for determining thermal 
conductivities are open to objection when used 


8 Eucken and Miicke, Zeits. f. physik. Chemie B18, 186 
(1932). 
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for polyatomic gases. When concentric cylinders 
are used instead of parallel plates, Eq. (d) has to 
be modified to 


—Q/2arNmD = C' (dT /dr) + Cyin(dTv/dr), 


where Q is the heat transfer/unit length of axis. 
This leads to integrals in the solution which have 
to be evaluated for definite values of the variables. 
The solution is not given here, however, since 
most methods with concentric cylinders require a 
correction for the temperature jump, particularly 
at the surface with the smaller radius of curva- 
ture. According to Gregory and Archer® the most 
satisfactory form of the correction for simple 
gases is to write AT’ /(zo+.S) for the temperature 
gradient in place of AT/zo. S is taken as pro- 
portional to 1/p. This correction is quite 
empirical and resembles the equation for F 
sufficiently to make a disentanglement of the 


® Gregory and Archer, Phil. Mag. 15, 301 (1933). 
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two effects somewhat problematical ; for evidence 
that this empirical correction alone is not 
adequate for CO: cf. M. Trautz.® 

Independent evidence of the connection be- 
tween thermal conductivities and the relaxation 
time for the various degrees of freedom can be 
obtained by considering the effect of small 
amounts of certain impurities. Thus the addition‘ 
of small amounts of He and H.O can increase 8 by 
several powers of 10; this should have the effect 
of increasing the thermal conductivity and 
decreasing the pressure coefficient.!° Such an 
effect has actually been observed" though here 
again the use of the vertical wire method makes it 
uncertain whether it is to be ascribed to con- 
vection effects, to a change of the accommodation 
coefficient in the presence of hydrogen, or to the 
effect discussed in this communication. 


10 Ubbelohde, Proc. Roy. Soc. A146, 271 (1934). 
1 M. Trautzand S. Weber, Ann. d. Physik 54, 498 (1917). 
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On the Definition of the Gibbs Potential 


J. N. BRONsTED, University of Copenhagen 
(Received February 1, 1935) 


N a short note On the Definition of the Gibbs 

Potential’ in which I commented on recent 
criticism by some authors as regards the in- 
troduction by Willard- Gibbs of the chemical 
potential into thermodynamics, I included the 
following remarks concerning an article by L. 
J. Gillespie and J. R. Coe, Jr. 2 “Quite recently 
thermodynamic relations of open systems have 
been deduced by L. J. Gillespie and J. R. Coe. 
They avoid intentionally the Gibbs potential as 
a basis for their deductions, on account of the 
fact that ‘Gibbs’ treatment has been held obscure 
by some authors.’ They say that Gibbs provides 
a physical interpretation (of the potential), but 


1 J. N. Brénsted, Det Kgl. Danske Videnskabernes Selsk. 
Mathem.-Fys. Medd. 12, No. 6, Copenhagen (1933). 
1933) Gillespie and J. R. Coe, Jr., J. Chem. Phys. 1, 103 


postpones it to a later part of his work.’’ As 
Professor Gillespie has kindly called to my 
attention that these remarks might be inter- 
preted as classing the two authors as critics of 
Gibbs, while the method adopted by them, 
avoiding the use of the potential, contained no 
such criticism, but was chosen for the sake of 
symmetry in the treatment, I am desirous to 
state that my remarks regarding the treatment 
of the two authors had no intention of such a 
classification. 

On the contrary I am aware that the physical 
interpretation of the chemical potential, which 
follows from a rigid analysis of the mathematical 
definition by Gibbs suggests itself also as a 
consequence of an analogous interpretation of 
the energy and entropy functions of open systems 
given by Gillespie and Coe. 
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The Resolution of Bond Eigenfunctions in Terms of a Linearly Independent Set 


HENRY EyRING AND HAROLD GERSHINOWITZ,! Frick Chemical Laboratory of Princeton University 
(Received December 27, 1934) 


By means of simple operators derived from the anti- 
symmetrical configuration eigenfunctions a method is 
devised for obtaining the coefficients in the expansion of a 
bond eigenfunction in terms of a linearly independent set. 
This method eliminates any necessity of drawing diagrams 
and removes the unique advantage of the uncrossed set. 
It is applied to the determination of the coefficients of 


eigenfunctions expanded in terms of the Rumer (uncrossed) 
set for as many as six bonds. For three and four bonds 
crossed sets are given for which the process of resolving 
eigenfunctions is simpler than for the Rumer set. A con- 
venient method of expansion in terms of the Rumer set is 
extended to include eigenfunctions of as many as six bonds 
of any multiplicity. 


I. INTRODUCTION AND GENERAL METHOD 


N order to take advantage of molecular 

symmetry in reducing the secular determi- 
nants of many atom problems it is necessary to 
be able to expand bond eigenfunctions in terms of 
a complete set. Until quite recently the calcula- 
tion of the coefficients in such an expansion 
offered considerable difficulty. A solution of this 
problem, complete from the theoretical point of 
view, has been given by Rumer, Teller and 
Weyl? who developed a geometric method of 
determining such a complete set and also show 
how by the laborious method of successive 
uncrossing of bonds, any eigenfunction can be 
expanded in terms of an uncrossed set. Their 
methods have been applied to the less compli- 
cated symmetrical molecules. As the number of 
bonds considered is increased the process rapidly 
reaches Herculean dimensions. Preceding papers 
from this laboratory have given rules for getting 
the coefficients in the expansions of eigen- 
functions containing as many as six bonds in 
terms of the Rumer set.* These methods, how- 
ever, enable one to use only one of the many 
possible independent sets. For systems of atoms 
having some other symmetry than that corre- 
sponding to rotations in a plane it is advan- 
tageous to use instead of the Rumer, that 
independent set which most nearly possesses the 
symmetry of the atomic configuration. We, 
therefore, considered it desirable to develop some 
method which would make more easy the 


1 Parker Fellow of Harvard University. 

2 (a) Rumer, Nach. Ges. Wiss. Gott. M. P. K. 1932, 332; 
(b) Rumer, Teller and Weyl, ibid. 1932, 499. 

5(a) Eyring and Sun, J. Chem. Phys. 2, 299 (1934); 
(b) Bear and Eyring, ibid. 3, 98 (1935). 


expansion of a bond eigenfunction in terms of any 
l.i.s. (linearly independent set). 

An analytical solution of this problem is 
suggested when one considers the function space 
of higher dimensions which underlies that of the 
bond eigenfunctions. A bond eigenfunction y; can 
be represented as a linear combination of the 
antisymmetric configuration eigenfunctions ¢; 
which have the same eigenvalue for S, the total 
z-component of spin. The ¢; are the determinants 
defined by* 


68(2), (1) 
Here N is the number of orbitals a, b, etc.; so that 
1/(N!)} is the factor which just normalizes the V! 
terms of (1) to the approximation of neglecting 
non-orthogonality of the orbitals. The quantities 
a and B are the two eigenfunctions for the S, fora 
single electron with the corresponding eigenvalues 
h/4x and —h/4r. The operation P permutes the 
electrons among the orbitals. Such an eigen- 
function can be represented symbolically in a 
very simple fashion, by indicating which spin 
eigenfunction is associated with each orbital. 
Thus for the singlet state of the four electron 


abcd 
problem we write ( a Bap for the eigenfunc- 


tion which has the eigenvalue S,=0, and which 
has 6’s in the 6 and d columns of its determinant. 

These ¢;’s are all linearly independent, whereas 
the y’s are not. Furthermore, the coefficients of 
the ¢;’s that combine to give a y; can be only +1 


4These are simply the antisymmetrical determinants 
used by (a) Slater, Phys. Rev. 34, 1293 (1929). For a more 
complete discussion see (b) Eyring and Kimball, J. Chem. 
He 1, 246 (1933); and (c) Bear and Eyring, reference 
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RESOLUTION OF BOND EIGENFUNCTIONS 


or 0, while in the expansion of a y; in terms of y's 
that are linearly independent, the coefficients are 
not thus restricted. Let us now suppose that we 
have an independent set Yi, and expand 
in terms of this set.5 


(2) 


r=1 


But we have also dix 


i=1 


(3a) 


l 
and (3b) 


i=l 
Therefore, substituting in (2) the expansions of 
y, and y, from (3a) and (3b) and equating the 
coefficients of 9; we get 


dia, 


r=1 


(4) 


It is convenient to think of the g’s (antisym- 
metric configuration eigenfunctions) as consti- 
tuting a set of mutually orthogonal unit vectors. 
The y’s (bond eigenfunctions) are then vectors in 
this space with one end at the origin and the 
other given by adding and subtracting vectorially 
an equal number of these orthogonal unit 
vectors. The coefficient 5;;, in (3a) is then just the 
scalar product ¢;:¥; (4a). This scalar product is 
not to be confused with the integral {y,*ydr, 
with which it agrees if { ¢;*9;d7 =6;;. The use of 
(4a) is of course not a limitation on our method of 
determining linear dependence. 

It is clear that bj, and 0;, can also be expressed 
as the simple operators D;= (0/0¢;) = , acting 
on y;, and y,, respectively, whence 


r=1 


(5) 


We evaluate Dw, for our independent set and 
substitute in (5), which process gives / equations 
(1>s) from any s of which we can solve for the s 
unknowns (the a,,) in terms of the Diy,. It is 
always possible to choose a set of s equations in 
such a way that one obtains linear equations of 
the form: 


5 In all that follows the subscript & refers to the specific 

nd function in which we are interested, and does not 
change, while r may be any one of the s linear independent 
bond functions and i any one of the / Slater functions. 
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The operator in the brackets is independent of k. 
We can therefore write 


(7) 
where (8) 


The operation A, on y;, of course amounts simply 
to taking a sum of scalar products. We have then 
a convenient method of obtaining the coefficients 
in the expansion of any bond eigenfunction. We 
need determine only once the operator A, for 
each y, and then the coefficient of that y, in the 
expansion of any y is given by A, operating on 
that y,. The simplest coefficients for any inde- 
pendent set will be found by picking the s 


. equations properly. In the applications that 


follow it will be shown how easily this may be 
done. 


II. APPLICATIONS TO SINGLET EIGENFUNCTIONS 
A. Three bonds 


Let us consider in detail the resolution of six 
electron bond eigenfunctions for the singlet state. 
For this case there are twenty ¢;’s and fifteen 
y,;’s, of which five of the latter are linearly 
independent. We shall arrange the ¢;,’s as follows. 
We shall write at the top of a column abcdef 
representing the six orbitals involved. Then we 
shall distribute three a’s and three §’s in all the 
twenty possible ways among these orbitals by 
arranging them as one would in constructing 
abcdef 


a dictionary. Thus ¢:= 


_ fabcdef _fabcdef , 


we arrange the ¢;’s in this manner we see that 
¢20 through ¢; are the reflections of ¢; through 
¢10, i.e., they have a’s where 8's were and vice 
versa. Thus the coefficients of goo through gi; in 
any expansion of a bond eigenfunction will be 
the negative of those of ¢; through ¢i9. Therefore 
in our further work we need only consider ¢; 
through gio. In the expansion of any bond 
eigenfunction y; the coefficients of the ¢;’s can be 
determined in the following manner. Corre- 
sponding to each ¢; we define an operator D; of 
the form [def] in which the letters simply list the 
orbitals in the corresponding y; which are 
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associated with 6’s. Operating on a bond eigen- 
function such as ¥4=ab cd ef, D; gives 0 if any two 
of the orbitals occur in any one bond; otherwise it 
gives +1 if the number of letters in it which are 
first letters of bonds is even and —1 if this 
number is odd. The following convention is to be 
understood for the ~’s. In each bond the letter 
coming first in the alphabet is written first and 
this eigenfunction we indicate by plus yx. 
Interchanging the order of the letters in one bond 
changes the sign of yx. 

Thus [def [bce ]Ws4=1; [bcf —1, etc. 
We consider, therefore, ten such operators, 
D,=([def], De=[cef], ---Dio=[bcd]. For our 
five linearly independent y’s let us first choose 
the Rumer set. 


vs=ab cd ef, 
vs=ab cf de. 


¥i=af bc de, 
bc ef, 
v3=af be cd, 


(9) 


Now any other eigenfunction y; can be written as 
(10) 


and we are interested in finding the aiz’s. 
Operating on y;, with the brackets defined above 
we obtain: 


diz =([def 

dex =([cef 

ds, =(edf 

dsx = (bef ¥x= 

dex = [bd f =G14— G24 O54, 
= 

dx = (bef 

do, = [bee 

ds ox = [bed = Asx, 


(10a) 


where 


dix (10b) 


From the five equations for diz, dex, daz, dox and 
dion we get the following solutions for the a’s. 


= dior. 


= — 4x, 
Q3k =dy, 


(11) 


For the expansion of any three-bond eigen- 
function y, in terms of the Rumer set we, 
therefore, have, 


—PoDitysDi 


lve. (12) 
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It is possible to pick an independent set in which 
not more than one operator is required to 
calculate each coefficient. This can be done by 
letting the ten D;’s operate on all fifteen y’s and 
then picking out five suitable y’s. When this is 
done we find that such a set is given by 
=ab cd ef, Yo2'=af be de, ce df, 
=ae be df, s'=af be cd. The coefficients 
are A;'=D,, A.’ =—D3, A;' Aj=—-D,, 
A;'=D,, where A; was defined in Eqs. (7) and 
(8) and D; was defined in (10b). The expansion 
of any y in terms of this set then becomes 


(13) 


Thus we see that it is a very simple matter to 
expand any bond eigenfunction in terms of any 
complete set. Even for the Rumer set it is no 
longer necessary to draw diagrams and repeated 
superposition patterns each time an eigenfunc- 
tion is expanded. Eqs. (12) or (13) need be 
obtained only once and then can be applied to 
any yx. 

Eq. (12) may be written in a simpler form by 
taking advantage of the symmetry properties of 
the Rumer set. The five eigenfunctions form only 
two patterns which can be represented by 
rotation (or permutation) groups of order two 
and three. Thus y is carried into — 4 by the odd 
permutation R=(abcdef) and into by 
repeating this permutation. The sign of an 
eigenfunction obtained by applying to a positive 
eigenfunction the permutation (abcdef) p times is 
(—)?. The same permutation changes y2 into 
—y3, into ¥;, and into As stated 
before the operation D; is simply ¢;-, so that the 
odd permutation R=(abcdef) applied p times 
gives =(—)?D;.° Furthermore, if the 
operator D; is one used to determine the coeffi- 
cient of y,, then (—)?R?D; is the operator 
associated with (—)?’R*y,. Thus we can rewrite 
(12) as 


6 This follows because any permutation applied toa ¢i 
gives some other Slater eigenfunction ¢;, and hence this 
permutation applied to the operator y;+ (which symbolizes 
the taking of a scalar product of g; with the eigenfunction 
that is operated upon) gives the operator (—)’—¢;:. 
which in our alternate notation is the operator (—)"Dj. If 
it should so happen that the permutation carries the given 
operator into one belonging to the set gu- to g20- inclusive 
(such as [aef], for example), one can either keep it or re- 
place it by the negative of its reflection from the set 
to gio (i.e. —[bed]). 
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2 
(14) 


This simple way of writing the expansion does 
not diminish the actual labor of carrying out the 
operation itself but will be found very convenient 
for the representation of eigenfunctions of more 
bonds where the number of those in the l.i.s. 
becomes quite large. Furthermore it does reduce 
the amount of labor needed to determine the 
operators since it is now necessary to find the 
operator for only one bond eigenfunction of each 
pattern. When the Rumer set is used it is found 
that those operators D; give the simplest 
coefficients ai, which have the largest number of 
6’s associated with neighboring orbitals (succes- 
sive letters). With these simplifications it becomes 
possible to determine quickly the expansions for 
eigenfunctions of many bonds. 

Another advantage of the notation employed 
in Eq. (14) is perceived when we proceed to 
determine the character of a representation. If 
we let some operator B act upon our independent 
set we get the character x of the corresponding 
transformation from the equation 


x=2A,Byr. (15) 
If to can be obtained from by the 
permutation P‘ and y,,; to ¥, by a permutation 
on we have 


x= TA BP + 


i=1 i=nt+1 


(16) 


In the particular cases for which B and P and B 
and Q commute it is seen that (16) can be much 
further simplified. 


B. Four bonds 


The Rumer set of four-bond eigenfunctions has 
fourteen members which fall into three groups of 
orders two, four and eight. These groups can be 
generated by the permutation R= (abcdefgh) 
operating on ¥,=ab cd ef gh; Y2=ad bc eh fg and 
¥3;=ab cd eh fg, respectively. The coefficients of 
these eigenfunctions are readily found to be 


= (Lbefg })¥i.=A x, 
=([edef 
a3, = ([bdef }+([cdef 


The eleven other coefficients are obtained from 
these by repeated application of the cyclic 
permutation R. We have then for the expansion 
of any four bond singlet eigenfunction in terms of 
the Rumer set, 


(17) 


FER (18) 


In the above set only four coefficients have one 
term, while eight have two terms and two have 
three terms. There exists a set for which eleven 
coefficients have one term, two have two terms, 
and one has three terms. It is the following: 


¥i1=ab cd eg fh, 
v2=ab ce dg fh, 
¥v3=ab cf dg eh, 
bd ef gh, 
¥5=ac be df gh, 
¥s=ad bg ce fh, 
vr =ad bg cf eh, 
vs =ac be dg fh, 
vo=ac af bg eh, 
vi0=af ce bg dh, 
vii =ad be cf gh, 
=ab ce df gh, 
¥i3 =ac bd eg fh, 
¥ia=ab cd ef gh, 


A, =[bdgh], 
A [begh], 
A;=[bfgh], 
A,=([cdfh], 
As=([cefh], 
Ao=[degh], 
A;=([dfgh], 
A = [cegh], 
Ao=[cfgh], 
A1o=([efgh], 
Ai =[defh], 
A13=[cdef ]—[cegh], 
A 


C. Five bonds 


There are 42 linearly independent five-bond 
eigenfunctions. The members of the Rumer set, 
however, form only six rotation groups, one of 
order two, two of order five and three of order 
ten. These six groups can be generated by the 
cyclic permutation from the eigenfunctions 


¥i=ab cd ef ghij 
¥2=ab cj di eh fg 
¥3=aj be di ef gh 
¥s=ab cj di ef gh 
vs=ab cj dg ef hi 
¥e=ab ej de fg hi 


(order 2), 
(order 5), 
(order 5), 
(order 10), 
(order 10), 
(order 10). 


The following are the coefficients of these eigen- 
functions from which the others may be gener- 


| 
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ated by the cyclic permutation R= (abcdefghij): 


ain = ([edefg]+ [bdefg ]—[bcdef 

+ [bedhi]+ (dfghi]+[efghi]) =A ivr, 
= [bedef Wx 
asx = (Lbcdef ]+ [bcdeg]+ 
4x = (Lbedef svi, 
asx = (Lbedef =A 
= (Lbcdef }+ (defgh]+[cefgh] 

— [bedgh]+ [befgh]}) =A 


D. Six bonds 


There are 132 linearly independent six-bond 
singlet eigenfunctions. The symmetry properties 
of the Rumer set reduce these to 14 rotation 
groups, one of order 2, one of order 4, three of 
order 6 and eight of order 12. We shall give in 
some detail the method of obtaining the coeffi- 
cients in this case as it throws some light on the 
way in which symmetry properties may be used 
to reduce enormously the labor of the calculation. 
Thus the operators for this set were obtained in 
several hours, whereas the method of uncrossing a 
single much crossed eigenfunction takes a time of 
an entirely different order of magnitude. Even to 
expand an eigenfunction by the simplified 
methods previously suggested’ would take longer 
than this operator method since in the former 
method to determine each coefficient it is 
necessary to draw a new diagram. 

Let us assume that twelve points a, b,c ---, 1, 
representing the orbitals, are arranged in order 
on the circumference of a circle (8). A bond 
between two of the orbitals is represented by the 
chord connecting the two points. Following a 
definition of Rumer, Teller and Weyl,’ we 
shall measure the length of the arc subtended by 
this chord, by the number of points on the arc, 
the two terminal points each counting one-half. 
We have already noted that the simplest coeffi- 
cients are given by those operators representing 
antisymmetric configuration eigenfunctions in 
which all the 8’s are on adjacent orbitals. We see 
at once that such an operator will give zero when 
acting on any bond eigenfunction except one of 
that symmetry group in which there are just two 
arcs of unit length. (Such an eigenfunction is 
ab cl dk ej fi gh.) The next simplest operator is 
that with five 6’s on adjacent orbitals. Such an 
operator will give zero when acting on all 
eigenfunctions except one of the symmetry group 


already mentioned and one of a group with just 
three arcs of length one. When one proceeds to 
choose the operators in such an orderly manner 
one finds that the solution of the fourteen 
simultaneous equations that are needed is made 
quite easy, since each new equation involved only 
the quantities that have occurred in previous 
equations. Thus, the fourteen eigenfunctions 
which follow are those of the Rumer set of 132 for 
which we require the corresponding operators. 


Order of rot. group 
¥i=ab cl dk ej fi gh 
v2=ab cl dk ej fg hi 
¥v3=ab cl dk eh fg ij 
v4=ab cl dk ef gj hi 
cl dk ef gh 1j 
ve =ab cl dg ef hk ij 
vz =ab cl di ef gh jk 
Ws=ab cl de fk gb ij 
9 =ab cl dg ef hi jk 
¥io=ab cd el fk gh ij 
¥u=ab cl de fi gh jk 
Yi2=ab cl de fg hi jk 
¥vis=ab cd ej fg hi kl 
Wis=ab cd ef gh aj kl 


The fourteen eigenfunctions given above are in 
the order of complexity of their operators. Thus 
the operator associated with y; has one term, that 
of ye two terms, etc. Furthermore, their sym- 
metry relations are such that each one contains 
the operator associated with the eigenfunctions 
that have gone before. This will be made clear by 
the following table. In this we give the operator 
coefficients of the eigenfunctions (A,), not 
explicitly in terms of the D,’s, but in terms of the 
A’s of the preceding eigenfunctions and whatever 
other D’s are necessary. In the actual expansion 
the labor is reduced if the coefficients are 
determined in this same order. This order can be 
followed with advantage in the cases of three, 
four and five-bond eigenfunctions. 

In order to get these relations in the simplest 
form it was necessary to use besides the operators 
for the fourteen eigenfunctions given above some 
of the others which are obtained by the repeated 
cyclic permutation. Thus in the expression for 
Ay the symbol R*‘A; means that one repeats four 
times the permutation R= (abcdefghijkl) on As. 
The expression R‘A; is the operator associated 
with the eigenfunction R43, which is ab ch dg ef 
al jk. 


— 


As 
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A; [bcdefg), 

Ag = [bedefh ]+A,, 

A; =[bedefi]+A2, 

A,=[bcdegh]+A2, 

As= [bcdegi ]— Ao, 

Ag=[bedehi]+As, 

A; =[bedegj ]— [bedefj ]+As, 

A 

Ag=[bedeij ]+A4+R'A;, 
Aic=[bcefhi]+As+ R8A7— 
Ai, 
Ai2=([bedghj ]+ [aefjkl]+ 
A 13 = [bcgijk—A itA R°A st % 
Aw = [bcfgjk]—A RA itA st+R?A 

+ R1As+ REA 5+ R8A 5 


III. THE EXPANSION OF BOND EIGENFUNCTIONS 
OF HIGHER MULTIPLICITIES 


The resolution of bond eigenfunctions of higher 
multiplicity could also be carried through by this 
same method of operators. However, for high 
multiplicities the number of linearly independent 
bond eigenfunctions becomes so large that the 
method begins to become unwieldy. For such 
cases one may use to advantage a combination of 
the geometric and analytic methods. The 1.i.s. for 
higher multiplicities may be obtained in the 
following way.” The orbitals a; to a, are 
arranged in order on a circle. Since we are con- 
cerned with higher multiplicities, not all of these 
orbitals are bound to others. A reference orbital is 
inserted on the circumference between a; and a». 
The independent set is that one composed of 
eigenfunctions in which no lines cross, when, 
beside the lines representing actual bonds, one 
adds additional straight lines connecting the 
unbound orbitals with the reference orbital. 

The first step in the resolution of bond eigen- 
functions of higher multiplicities is to disregard 
all orbitals that are not connected by bonds and 
to expand the resulting eigenfunction in terms of 
an uncrossed set by the methods of the preceding 
section for the singlet state. These resulting 
uncrossed eigenfunctions will not in general be 
members of the I|.i.s. We must, then, expand each 
of them. 

It has been shown’ that there is a vector-like 
relationship between eigenfunctions which are 
exactly the same except for a single bond. If 
these differing bonds when superposed on the 
same circle form a closed polygon any one of the 


7 Eyring, Frost and Turkevich, J. Chem. Phys. 1, 777 
(1933), 


eigenfunctions can be represented as the vector 
sum of the rest. This principle, together with the 
definition of a l.i.s. given above, enables us to 
expand our uncrossed eigenfunction in terms of 
the L.i.s. 

The method of expanding an eigenfunction of 
any number of bonds and of any multiplicity is 
now summarized. After one has applied to a 
crossed eigenfunction the method for obtaining 
an uncrossed set, as previously described, each 
one of the set is separately expanded. First we 
introduce some definitions. We speak of the bond 
which has one of its ends farthest from our 
reference orbital as By the one next to it as By, 
etc. By the arc subtended by B, is meant the arc 
which does not contain the reference orbital. We 
next define a quantity which we will call the 
weight of a bond. The weight of By is equal to the 
length of its arc mo (defined previously). The 
weight of B; is its length minus 27. The product of 
the weights for all bonds of an eigenfunction y 
gives the number of eigenfunctions of the L.iss., 
which will occur in the expansion of y. If for 
every eigenfunction the orbital in a bond with the 
lower ordinal number is written first all these 
eigenfunctions will occur with a coefficient +1. 
This is the convention of sign which has been 
employed throughout the paper. | 

To resolve the eigenfunction y one first 
expresses it as the sum of eigenfunctions which 
differ from y only in having By replaced in turn 
by each possible bond B,’ lying along the arc 
subtended by Bo. Each of the resulting new 
functions is in turn replaced by a_ sum. 
This sum is obtained by replacing B,, suc- 
cessively, by all the bonds which add vectorially 
to give B, and which follow along the arc of B, 
at all points except at the two points corre- 
sponding to the extremities of the By’ for this 
eigenfunction. These points are omitted. This 
same process is continued for all the bonds up to 
B,_; (there being assumed to be m bonds). In 
expanding a particular eigenfunction by suc- 
cessive replacements of the bond B; we proceed as 
for B; except that we now omit the 27 points on 
the arc of B; which coincide with the two 
extremities of the 7 bonds By’ B,’---B’;_;. Thus 
eigenfunctions for higher multiplicities, for any 
number of bonds n, are readily resolved in terms 
of the linear independent set. 
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The Quantum-Mechanical Treatment of Molecules by the Method of Spin Valence* 
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The method of spin valence, developed by Born, Heitler, 
Rumer and Wey] for the treatment of molecular systems, 
becomes completely equivalent to the Slater method if the 
bonding forces are considered to be exerted between in- 
dividual orbitals, and not between entire atoms in definite 
spectroscopic states. This change in the point of view 
eliminates several difficulties which were formerly present 
in the method (e.g., the obscuring of the directional proper- 
ties of the valence bonds). The secular equation, for which 
a simplified derivation is given, is of quite a different form 


from the Slater equation, and the problem of its solution is 
considered in some detail. In addition to the rigorous 
methods, several approximate methods are given, of which 
two (applicable to unsaturated and aromatic hydrocarbons) 
do not require that the secular equation be set up at all. 
The method of spin valence is found to be more convenient 
than the Slater method for the treatment of systems of high 
multiplicity. Finally, the extension to cases of orbital, as 
well as spin degeneracy is briefly outlined. 


INTRODUCTION 


HE method outlined by Slater! for the 
quantum-mechanical treatment of molec- 

ular structure is of such general applicability that 
certain other procedures which have been devised 
and which lead to the same final result are apt to 
be completely neglected. (I am not concerned 
here with the method? associated chiefly with the 
names of Hund and of Mulliken, since this starts 
from an essentially different basis and is not 
strictly comparable with Slater’s treatment.) 
- Van Vleck and Serber,* to be sure, have recently 
pointed out that the Dirac vector model is quite 
useful for many purposes, but very little atten- 
tion seems to have ever been paid to the so-called 
method of spin valence developed by Born, 
Heitler, Rumer and Weyl.‘ One reason for this 
neglect is probably to be found in the fact that, 
in its usual formulation, the method treats only 
the interactions between entire atoms in definite 
spectroscopic states, and not those between the 
individual electrons. This makes it almost neces- 
sary to assume each atom to be in an S state, and 
obscures the directional properties of the valence 


* No. 455. 

1]. C. Slater, Phys. Rev. 38, 1109 (1931). 

2 F, Hund, Zeits. f. Physik 73, 1 (1931); 73, 565 (1932). 
R. S. Mulliken, J. Chem. Phys. 1, 492 (1933) and further 
references given there. Cf. also E. Hiickel, Zeits. f. Physik 
72, 310 (1931); 76, 628 (1932); 83, 632 (1933). 

3 J. H. Van Vleck, Phys. Rev. 45, 405 (1934). R. Serber, 
Phys. Rev. 45, 461 (1934); J. Chem. Phys. 2, 697 (1934). 

4M. Born, Zeits. f. Physik 64, 729 (1930); Ergebnisse 
d. Exakt. Naturwissenschaften, X 387 (1931). W. Heitler 
and G. Rumer, Zeits. f. Physik 68, 12 (1931); H. Weyl, 
Nach. Ges. Wiss. Gott., M. P. Klasse, 1930, 285; 1931, 
33; G. Rumer, ibid. 1932, 337. 


bonds to such an extent as to invalidate many of 
the conclusions reached. The incorrect result 
obtained by Markov,' for example, that benzene 
should be less stable than three acetylene mole- 
cules is apparently to be explained on this basis. 
This difficulty, however, is not at all essential, 
and, as we shall see immediately, can easily be 
avoided. 

In the following development, a method of 
attack will be adopted which differs from that 
usually employed in that it assumes the Slater 
valence-bond functions at the outset. This 
eliminates at once the above-mentioned difficulty, 
and simplifies the subsequent treatment by 
rendering unnecessary any recourse to the 
unfamiliar theory of the invariants of binary 
forms. With this alteration, the method of spin 
valence becomes completely equivalent to the 
ordinary Slater treatment, although it remains 
quite different in form. It always leads to the 
same answer as the latter in any definite problem, 
but it frequently allows the calculation to be 
made with the expenditure of considerably less 
time and labor. Even in the cases in which it 
offers no special advantages, it is still useful in 
that it furnishes an independent check upon the 
results obtained in other ways. 


OUTLINE OF THE METHOD 


For the present, the discussion will be limited 
to problems of spin degeneracy only, and to 
systems in singlet states. If we have 2N orbitals, 
a, b, -++y, 2, each occupied by one electron, we 


5M. Markov, J. Chem. Phys. 1, 784 (1933). 
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can set up a Slater function for each possible way 
of drawing valence bonds. For bonds between a 
and 6, c and d, ---y and z, for example, we have® 


1)? 


X {a(1)e(1)b(2)B(2) - -9(2N—1) 


Xa(2N—1)2(2N)8(2N)}, (1) 


where P; is one of the (2N)! operations of 
permuting the electrons among the orbital and 
spin functions, and R; is one of the 2% operations 
of interchanging the spins associated with the 
orbitals involved in one or more bonds. P; is an 
even, or an odd, permutation when 7 is even, or 
odd, respectively; and R; corresponds to an even, 
or to an odd, number of interchanges when b is 
even, or odd, respectively. In order to avoid 
confusion in regard to the sign of the function, we 
shall adopt the following convention: Let each 
valence bond in the structure be represented by 
an arrow. Then the corresponding function is 


X 


where 


oa t™ if the orbital 7 is at the tail of an arrow 
"| 8, if the orbital 7 is at the head of an arrow. 


In that term of the double summation for which 
R,, and P; are both equal to the identity operator 
I, the orbitals follow in alphabetical, and the 
electrons in numerical order. Reversing an arrow 
in the structure then changes the sign of the 
corresponding function. In the structure repre- 
sented by Eq. (1), for example, the arrows point 
from a to b, ---, from y to z. 

If the summation over k is actually performed, 
Eq. (1) becomes 


1)*P;{a(1)b(2)--- 


X2(2N)[a(1)6(2) 


° The first explicit statement of this rule (in a slightly 
different, but equivalent form) seems to have been made 
by G. E. Kimball and H. Eyring, J. Am. Chem. Soc. 54, 
3876 (1932). 


[a(2N — 1)8(2N) —B(2N—1)a(2N) }} 


(—1)i(Pju)(P; 
)(Pju)(Pi¢) 


where, for convenience, we have set 


u=a(1)b(2)---2(2N) 
and 


y=[a(1)6(2) —B(1)a(2)]--- 
X [a(2N — 1)8(2N) —B(2N—1)a(2N)]. 


From this original structure, any other that 
may be desired can be obtained by redirecting 
the arrows in a suitable manner. Such a redirec- 
tion can be represented by a permutation, 7;, 
acting upon the orbitals. If, for example, 7’; is the 
simple interchange (by), the corresponding struc- 
ture is that in which the arrow, which formerly 
had its head (or tail) at b, now has its head (or 
tail) at y, and vice versa. That is, 


becomes 
y 


Any other bonds (such as between ¢ and d, etc.) 

remain unchanged. It is readily seen that the 

function which represents this new structure is 

(Pye). 

i 
(2%(2N)!)¥ 


a—b 


Z—y 


(4) 


The factor (—1) ensures that the direction of 
each arrow in the vector-bond structure will 
remain unchanged in the course of the per- 
mutation. 

We can put Eq. (4) into a more convenient 
form by noting that 7; can be replaced by P;"", 
where P; is the same permutation of the elec- 
trons, 1, 2, ---2N, as T; is of the orbitals, a, b, 

-+g. This gives 


i+i(P; 
(Pig). (5) 


If now we denote P;P ;— by P;, we have, since P; 
and P;~ are even or odd simultaneously, 


1 
*(Pyu)(P 


or, dropping the primes and setting Pig= ¢; 
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1 
Pu) (Pigs). (6) 


(2%(2N)!)8 i 

There are (2N)! permutations 7; and hence 
(2N)! functions ¥; which can be obtained in this 
way, but of these only (2N)!/[N!(N+1)!]= pw 
can be independent, since there are only py 
singlet states. It is obvious, for example, that if 
T; is the interchange (ab) in Eq. (4), the effect is 
merely to change the sign of the function, or to 
reverse the arrow between a and 3 in the corre- 
sponding structure. The general case can be 
described in the following manner. Arrange the 
orbitals formally in a circle. (This need have no 
relation to the actual arrangement in space.) 
Now draw all possible structures in which each 
bond lies wholly within the circle, and in which 
no. two bonds intersect. For this purpose struc- 
tures which differ only in the direction of one or 
more of the arrows are not to be considered as 
distinct. There will be exactly py such structures, 
and the corresponding functions will form a 
complete and independent ‘‘canonical”’ set. 


Structures which contain crossed bonds can: 


always be resolved into linear combinations of 
the canonical structures by repeated application 
of the rule’ 


X= = - 


or by the procedure given by Eyring and Sun.® 

In general no single function WV; will be a 
satisfactory representation of the molecule, and it 
is necessary to use a linear combination 


x=L a; 


i=1 


(7) 


in which the a;’s are constants which are to be so 
chosen that x is the best possible approximation. 


Instead of using the variational method, which. 


leads to Slater’s secular equation, we shall 
evaluate the coefficients by the ordinary first 
order perturbation theory. Thus, we set 


(H-W)x= a(H-W)%i=0. (8) 


If we should multiply Eq. (8) by ¥, (J=1, 2, 
-++by) and integrate over all coordinates, we 
7G. Rumer, Nach. Ges. Wiss. Gétt., M. P. Klasse, 1932, 


337. 
8 H. Eyring and C. E. Sun, J. Chem. Phys. 2, 299 (1934). 


would again end up with Slater’s secular equa- 
tion. In the present treatment, however, we 
multiply by u, and integrate only over the spatial 
coordinates. In this way we obtain: 


X (u| 
or if we set 
G;=(—1)i(u/H— W/P ju) 


Pn 
Pjei=d gi, 


l=1 


and 


Py PN 


YX a yi =0. 


i=l l=1 


we have 


The ¢,’s are linearly independent, and conse- 
quently the coefficient of each must vanish 
separately. Thus 


2. (i) Ja: = 0, 


i=1 7 


(J=1,2,---pw). (9) 


This set of simultaneous homogeneous linear 
equations for the a;’s has a solution if and only if 
the determinant of the coefficients vanishes: 


UGibu LG by 
7 7 
=0. 
> Gjboy ty 
7 


(10) 


7 


This secular equation is quite different in form 
from the Slater secular equation which it re- 
places, and its matrix elements are obtained in 
quite a different way, but it clearly must lead to 
the same values of the energy W, and of the 
coefficients a;. It is particularly to be noticed 
that the determinant is not in general sym- 
metrical about the main diagonal, since )°G,b;)” 


is not necessarily equal to )}Gjbi;5. 
7 


AN EXAMPLE 


In order to illustrate the method, let us con- 
sider benzene, which has already been treated by 
Hiickel,® and by Pauling and Wheland.!° We 
shall take into account only the interactions 
between the six [p], orbitals (using Hiickel’s 

* E. Hiickel, Zeits. f. Physik 70, 204 (1931). 


10L. Pauling and G. W. Wheland, J. Chem. Phys. 1, 
362 (1933). 
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TABLE I. 


(bc) (cd) (de) 


(fa) z 


—A —D+A —A 


A(Q—W)+(C+D+E)a 


—C+B —B —E+B 


—-D+B B(Q—W)+(C+D+E)a 


-c | | -A+C 


—-A+C C(Q—W-—2a)+2(A+B)a 


—-A+D —D —-A+D 


—B+D -D D(Q—W—2a)+2(A+B)a 


—-A+E | —B+E —E 


—-B+E | -—A+E E(Q—W-—2a)+2(A+B)a 


nomenclature) and we shall make the usual 
assumption that G;=Q—W if T; is the identical 
permutation J, G;j=—a if 7; is a single inter- 
change between adjacent orbitals, and G;=0 
otherwise. There are five canonical structures, 


a a 


a 

b f fob 

d 


d d 
D E 


The actual calculation is given in Table I. 
At the top are listed the permutations considered 
(that is, those for which G;#0), and on the left 
are listed the independent spin functions. For 
convenience, ga is written simply A, and so on. 
The Slater function which represents the struc- 
ture A, for example, is of course 


= (P; P; 
Va 


This gives W=Q-—2a, Q—2a, Q—(13!+1)a, 
Q+(13!—1)a, and consequently the roots are 
exactly the same as those found by Pauling and 
Wheland from the Slater secular equation, and by 
Hiickel by a much more complicated procedure. 


In each square of the body of the table is given 
P;; corresponding to the permutation 7; above 
and to the spin function g; at the left. As an 
illustration, the permutation 7;= (ab) transforms 
the structure A into 


a 
f b 
d d 
Thus, 7;A=A-—E, and, consequently, Pj¢a 


= g4—¢e. In the column at the extreme right of 
the table are the values of 


E 
i 


j l=A 


obtained by summing the corresponding row. 
From these entries the secular equation can be 
written down at once: 


Q-—-W-2a 


The eigenfunction corresponding to the lowest 
root, W=Q+(13!—1)a, is before normalization 
Wat+Vp4+0.4341(Vce+Wp+WVez). This is again the 
same result as was obtained by Pauling and 
Wheland. 
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T; 
A | A | -E+A —C+A 
clc|-3(| | 
| D | —B+D = 
a a 
f b f b f 
t t 
d d 
A B 
a a a 
0 Q-W a a a 
2a 2a Q-W-2a 0 0 =0. (11) 
2a 2a 0 Q-—W-2a 0 
2a 0 


234 


UsrE oF SYMMETRY OF THE MOLECULE IN SIMPLI- 
FYING THE TREATMENT 


The calculation as carried out above is of about 
the same difficulty as that by the Slater treat- 
ment, but it can be greatly simplified if the 
energy of only the lowest state is desired. This 
follows from the fact that the degree of the 
secular equation can be reduced by taking 
advantage of the symmetry of the molecule." It 
is obvious, for example, that the structures A 
and B must enter the complete eigenfunction for 
the ground state with the same coefficient, and 
that C, D and E must do likewise. (In order for 
this statement to hold, it is necessary that each 
arrow in each structure point from an odd 
numbered to an even numbered orbital, the 
orbitals being numbered in sequence around the 
circle in either direction. The starting point for 
the enumeration is arbitrary, but must be the 
same for all structures. The purpose of this 
convention is toensure that the integral (VW ;W;) >0 
for all 7 and 7.) There are thus only two inde- 
pendent coefficients, and in consequence the 
secular equation can be made a quadratic. It is 
not necessary, however, to set up the original 
fifth degree equation at all, since it is possible to 
obtain the quadratic in one step. The procedure 
in the general case is as follows: Let there be 
several different types of structure, of which 
characteristic representatives are A, B, ---F, 
respectively, and let there be ga structures of 
type A, all with coefficient a4, gz structures of 
type B, all with coefficient ag, and so on. Of 
course g4+ge+-°:+++gr= py, the total number of 
independent spin functions. Now, using only the 
structures A, B, ---F as a basis, construct a 


11 This is a special case of quite general group theoretical 
considerations. Cf. H. Eyring, A. A. Frost and J. Turke- 
vich, J. Chem. Phys. 1, 777 (1933); A. E. Stearn, C. H. 
Lindsley and H. Eyring, ibid. 2, 410 (1934). 
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table similar to Table I. If the carrying out of a 
permutation upon a given member of the set 
A, B, ---F leads to a structure not a member of 
the set, this latter is represented in the table by 
that member of the set which belongs to the 
same type. The secular equation is then written 
down from the table in the same manner as 
before. If the coefficients a4, ap, -+*apr are de- 
sired, they can be obtained from the equations 


P 
(Mrs—srsW)argr=0, (S=A, B, +++F), (12) 


R=A 


where Mrs—irsW is the RS element of the 
secular equation. These equations take the place 
of (9), from which they can be readily derived by 
equating all the a;’s that should be equal, and by 
then adding together the equations for which / 
refers to structures of the same type. 

As an application of this procedure, we shall 
consider cyclooctatetrane, which has already 
been treated by Penney” by the Slater method. 
There are 8 electrons and hence 14 canonical 
structures, of which 2 are of type A below, 8 of 
type B, and 4 of type C. The remaining structures 
of each type are obtained from the characteristic 
example by rotation about the eightfold axis. 


fre fre 
B Cc 


The results of the different permutations are 
given in Table II. It is understood of course that 
the entry corresponding to T;= (de) and 
for example, does not mean that the result of the 


12, W. G. Penney, Proc. Roy. Soc. A146, 223 (1934). 


TABLE II. 


(fg) 


A(Q—W)+4Ba 


—A+B B(Q—W)+(24+C)a 


x 

¢ 
A 
B 
Cc 


—B+C C(Q—W—4a) 


+(2A +4B)a 


st 
T 
T 
tk 
at 
ec 
Se 
de 
a] 
ql 
cc 
th 
ab a b a<b or 
h é h lit 
g d 
sh 
A al 
pl 
is 
sy 
cc 
| I | (ab) (be) (cd) (de) (ef) 
Ci =-€ | -84C1 
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permutation (de) upon C is actually equal to 
—B+C, but merely that it is equal to —(a 
structure of type B)+(a structure of type C). 
The secular equation is now 


—12(Q— W)a?+40a' =0. 


This gives VW= Q+3.3022a for the lowest root, in 
complete agreement with the value found by 
Penney. The work involved here is much less 
than by the Slater treatment, which necessitates 
actually setting up the entire fourteenth degree 
equation, and then reducing it to a cubic. 


APPROXIMATE SOLUTION OF THE 
SECULAR EQUATION 


First method 


For still more complicated molecules, the 
secular equation becomes of unmanageably high 
degree, even after all the reduction arising from 
symmetry has been taken into account. In such 
cases it is desirable to have some method of 
approximation. In the Slater treatment this is 
quite easy, since one has only to equate the 
coefficients of all the first excited structures, 
those of all the second excited structures, and so 
on. One can furthermore entirely neglect the very 
highly excited structures since these have but 
little effect upon the total energy. Sherman has 
shown!® by a detailed discussion of naphthalene 
that this procedure results in only a small error, 
and is a very good approximation indeed. In the 
present treatment, however, the secular equation 
is of quite a different form. In particular, it is not 
symmetrical about the main diagonal, and, as a 
consequence, the sign of the error which is 


‘J. Sherman, J. Chem. Phys. 2, 488 (1934). 
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introduced by the various approximations cannot 
be known beforehand with any certainty. 
Furthermore, the precise significance of the 
energy value obtained when the equation is not 
solved rigorously is no longer clear. Conse- 
quently, we cannot be sure, without further 
discussion, that the same approximate procedure 
as was used with the Slater equation will be 
legitimate here, but we can apply it, nevertheless, 
to several problems, and see how well it works in 
practice. This will give only an empirical check 
upon the method; but a rigorous theoretical 
treatment would be quite difficult. 

As a first example, let us consider cyclo- 
decapentane. With 10 electrons there are 42 
canonical structures, which belong to the 6 types 


gp=10 ge=10 


The secular equation is found in the usual manner 
to be: 


| 
a 
t 
of 
iS 2a 4a Q-—W-—4a 
/ a—>b a—b b 
N 
j j<—— j 
Il 7 
y d i do > 
1 h e h e h e 
f XY \ 
g f gf g f 
c A B Cc 
3. ga=2 ge=10 
a 
jt j 
7 
6 
eof 
: D E F 
gr=5 
QO-W 5a 0 0 0 0 
2a Q-W-2ea 2a 0 2a 0 
- 0 4a Q-—W-2a 2a 0 0 
(13) 
2a 0 2a Q—-W-—4a 2a a 
0 Sa 0 2a 0 
0 4a 0 4a 0 
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This equation, being of the sixth degree, is still 
not too complicated to be treated rigorously. On 
multiplying out the determinant, we obtain 


(Q— W)*—18(Q— W)*a+84(Q— W)4a? 
+104(Q— 1520(Q— W)?at 
+3104(Q— W)a®—1600a°=0 


of which the lowest root is W=Q+4.031a. The 
correct value will serve as a useful check upon the 
approximate treatment, which we now apply. 

Before we can do this, however, we need to 
observe that while the secular equation enters 
the calculation as the condition of compatibility 
of the set of Eqs. (12) for the argr’s, it can also 
be considered as the condition of compatibility of 
a different set 


(Mrs— Wérs)cs=0 (R=A,B,-:-F). (14) 


S=A 


In the Slater treatment, the constants cs are just 


2(Q-—W) 10a 
20a 15(Q—W) +10 

20a 70a 

0 20a 


when we group together the first excited struc- 
tures B and C, and the second excited structures 
D and E. On multiplying out the determinant, we 
find 


100[3(Q— W)*+ 22(Q— W)8a—27(Q— W)?a? 
— 348(Q— W)a3—380a*]=0 


or W=Q+4.022a, for the lowest root. This value 
is in error by only 0.009a, and consequently is 
very good. If now we introduce the further 
simplification of neglecting the third excited 
structure F, we have merely to solve the cubic 
equation obtained by dropping the fourth row 
and column from (15). This leads to an energy of 
W=Q+3.954a, in error by 0.077a, which is still 
fairly satisfactory. If, however, the second 
excited structures, D and E, are neglected as 
well, the energy obtained from the resulting 
quadratic equation is only W= Q+2.937a; and if 
all excited structures are neglected, the energy is 
simply W=Q. Thus the neglect of the highly 


20(0— W) —40a 10a 


equal to the corresponding as’s, since the secular 
equation is symmetrical, but in the present 
treatment, this is not necessarily true. The c’s 
here are in general quite complicated functions of 
the a’s, and their precise significance is not at all 
clear. Consequently we have no intuitive reason 
(as we do in the case of the a’s) to believe that 
the c’s for structures of the same degree of 
excitation should be approximately equal to each 
other. We assume, however, that they are, and 
check the legitimacy of the assumption by 
comparing the result thus obtained with the 
rigorous solution. We are justified in doing so, 
since the present approximate treatment is 
almost entirely empirical and exploratory in 
nature. 

The procedure now is first to multiply each 
row of the secular equation by the corresponding 
Zr value," and then to add the rows and columns 
which refer to structures whose coefficients are to 
be equated. Thus the Eq. (13) becomes 


0 0 
30a 0 


20a 5(Q—W) —30a 


excited structures is not as satisfactory here as in 
the Slater treatment, and must be employed with 
considerable caution. 

These same general conclusions are reached in 
a consideration of naphthalene. Here there are 


again 10 electrons and 42 canonical structures, . 


but the symmetry of the molecule allows the 
reduction of the secular equation to be carried 
only as far as the sixteenth degree. Sherman’ has 
given a complete treatment of this problem by 
the Slater method, and has shown that the 
rigorous solution is W=Q+4.0400a. We shall 
not solve the present equation rigorously,” but 
shall reduce it to the fifth degree by equating the 
coefficients of the first excited structures, those of 
the second, and those of the third. This gives 


4 The purpose of this step is to ensure that the subse- 
quent adding of rows will equate the coefficients az, and 
not the products argr. 

16 Dr. J. Sherman and the author, however, have proved 
that it is completely equivalent to the Slater equation. 
The proof rests upon the fact that it can be made identical 
with the latter by the use of the familiar transformations 


d 
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(The same energy would be obtained if this 
equation were reduced still further to the fourth 
degree, by grouping the two different types of 
unexcited structures; i.e., by adding the first two 
rows and columns.) The lowest root is found to be 
W=Q+4.0241a, in error by only 0.0159a. This is 
somewhat better than the value of W=Q 
+4.0175a which Sherman obtained by applying 
the same approximate treatment to his Slater- 
type equation. The neglect of the third excited 
structures leads to an energy of W=Q+3.7231a, 
which is to be compared with the value of 
W=Q+4.0172a, obtained by Sherman with the 
same assumption. Again we see that this latter 
approximation is not very satisfactory in the 
present method. 


Second method 


It is readily seen by comparing (12) and (14) 
that the equation 


(16) 


is rigorously true for any arbitrary choice of the 
constants cs if the a’s are correct, and for any 
arbitrary choice of the constants ae if the c’s are 
correct. Consequently, it must be very nearly 
true if both the a’s and the c’s are approximately 
correct. Rough values of the various coefficients 
can be obtained for this purpose most easily by 
the method of successive approximations. That 
is, we first arbitrarily insert in the Eqs. (12)—or 
(14)—any values of the a’s—or of the c’s—that 
may seem reasonable. In this way, we obtain a 
set of linear equations for W, the roots of which 
will not be mutually consistent unless the 
assumed values of the coefficients are accidentally 
correct. We then proceed systematically to 
improve the consistency by varying the coeffi- 
cients one or two at a time. With a little practice, 


of determinants (i.e., adding columns, etc.). This gives a 
Practically complete check upon both equations. 
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Q-—W-a 0 6a 0 
0 2(0—W) —2a 12a 0 
14a 20a 16(Q— W)—10a 38a 2a =0. 
8a 12a 72a 19(0— W) —47a 12a 
4a 8a —4a 20a 4(Q—W)—20a 


this becomes quite a simple operation, and no 
great labor is required to obtain fairly satis- 
factory values. For example, in the case of 
cyclodecapentane, we find very easily that 


p21:1/3:1/4:1/15:1/10:1/80, 
and 
2 1:4/5:4/5:18/25:18/25:3/5. 


When these values are inserted into Eq. (16), 
the resulting energy, W=Q+4.037a, is in error 
by only 0.006a. 


FURTHER APPROXIMATE METHODS APPLICABLE 
TO UNSATURATED AND AROMATIC 
HYDROCARBONS 


If no great accuracy is desired, this last 
procedure can be made into a very simple one by 
the use of merely estimated coefficients. A fairly 
satisfactory empirical rule for doing this with the 
hydrocarbon molecules is as follows: Let a=c=1 
for all unexcited structures; a= 3/10 and c=5/6 
for all first excited structures; a= (3/10)? and 
c= (5/6)? for all second excited structures; and 
so on.* The sort of results that can be obtained by 
this method is indicated in Table III (third 
column from the left). 

A still rougher approximation consists in 
letting all the c’s be equal to unity. Then Eq. (16) 
becomes 


or 
¥ argeMrs/>d anger. (17) 
R S R 
Now, if we examine the various secular equations 
which have been set up above, we notice a very 
striking regularity, namely, that in all cases 


* Note added in proof: Further investigation has shown 
that this rule is greatly oversimplified. The author expects 
to discuss the problem in further detail in al ater paper. 
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TABLE III. Approximate energies of unsaturated and aro- 
matic hydrocarbons. 


Energy 
Compound Correct*® Approximate 
II4 


Benzene Q+2.61la Q04+2.5a 
Cyclooctatetrane 3.30a 34a 3.3a 
Naphthalene 
Butadiene 
Hexatriene 
Octatetrane ~ 
Phenylethylene 
Diphenylethylene 
Triphenylethylene 
Tetraphenylethylene 
Biphenyl 
Diphenylbenzene 
Triphenylbenzene 


Sa 


®L. Pauling and J. Sherman, J. Chém. Phys. 1, 679 (1933). 

b These values were obtained by the Slater method, with only un- 
excited and first excited structures taken into consideration. Conse- 
—_ the coefficients of a appearing here are lower limits of the true 
values. 

¢ Values in this column were obtained from Eq. (16) with the use of 
estimated coefficients: a =(0.3)i and c =(5/6)i for all ith excited struc- 
tures. 

4 Values in this column were obtained from Eq. (19) with a's esti- 
mated as in column I. 

¢ Values in this column were obtained from Eq. (21). 


Mrs=Q-+ne, (18) 


where 7g is the number of strong bonds between 
adjacent orbitals in the structure R. That is, 
nr=N for the unexcited structures, 
for the first excited structures, and so on. This 
rule is easily shown to be quite general, as long 
as the only interchanges to be considered are 
those between odd numbered and even numbered 
orbitals, the enumeration being carried out as 
described above. (page 234.) (In the type of 
problem under consideration, this is equivalent to 
the restriction that only interchanges between 
adjacent orbitals be considered, and that the 
molecule contain no odd-membered rings.) Thus, 
putting Eq. (18) into (17), we find 


W= Q+d (19) 


This relation makes it possible to calculate the 
approximate energy of the molecule without 
setting up the secular equation at alJ. One needs 
only to know the factors ge and the approximate 
values of the coefficients ar. Some results ob- 
tained in this way—with the same assumptions as 
before in regard to values of the a’s—are given in 
Table III (second column from the right). It will 
be seen that these are really quite good. The 
value of this method—and of the immediately 
preceding—lies in the fact that the calculated 


energy is not very sensitive to small changes in 
the coefficients, as long as these are of reasonable 
magnitudes. 

Since Eq. (19) is rigorously true when the 
correct values of the a’s are inserted, it gives a 
very valuable check upon the results obtained by 
other methods. For example, when the coeffi- 
cients given by Sherman! (column 1a in Table 
III, page 491) for the 16 different types of 
structure of naphthalene are used, the resulting 
energy is W=Q+4.0403a, in almost complete 
agreement with his value of W=Q+4.0400a. 
Consequently one can be quite certain that there 
is no large error in his work. When the degree of 
the secular equation is reduced by equating 
coefficients arbitrarily, Eq. (19) is no longer 
necessarily exactly true. Thus, with Sherman's 
fifth degree coefficients (column 1b in Table III), 
the energy becomes W=Q+4.0545a, which is 
different from the solution of the corresponding 
equation—W= Q0+4.0175a. 

There is one final approximate method which 
is useful for order-of-magnitude calculations. 
Like the last, it does not require that the 
secular equation be set up at all, and, in addition, 
it does not require any knowledge of the factors 
gr. If all unexcited structures are given the 
coefficient c=1, and all first excited structures 
the coefficient c= ¢, then one can easily show that 
any equation of the set (14), for which R refers to 
an unexcited structure, can be expressed in the 
form 


(20) 


Here B is defined as the number of interchanges 
in the molecule for which the corresponding 
exchange integral is equal to a—i.e., 6 for 
benzene, 11 for naphthalene, and so on. (The 
derivation of (20) assumes that all unexcited 
structures can be canonical simultaneously, and 
also that there are no ‘“‘phantom orbitals,”’ as in 
the case of the free radicals.'®) The extreme 
variation of @ seems to be approximately from 
0.75 to 0.9 (this would give an upper and a lower 
limit’to the energy), although it is possible that 
this limitation may not be valid in molecules 
widely different from those studied so far. If a 


rounded-off average value of €=5/6 is used, Eq. 


6 Cf. L. Pauling and G. W. Wheland, J. Chem. Phys. 1, 
362 (1933). 


e- > = 


| 
i 
0 
c 
t 
v 
7 il 
b 


SPIN VALENCE 239 


(20) becomes 

W=Q+ 6(7N—B). (21) 
The results of a number of calculations by this 
method are given in the last column of Table IIT. 
It will be seen that these are considerably better 
than one would have any right to expect. 


STATES OF HIGHER MULTIPLICITY 


The method of spin valence is especially 
advantageous in the treatment of systems of high 
multiplicity. The Slater method, of course, can 
also be extended to these cases, but, as far as 
the author is aware, no very convenient procedure 
for doing so has yet been described. The device of 
formally reducing the system to a singlet by the 
introduction of ‘“‘phantom orbitals’’!? becomes 
very cumbersome as the multiplicity increases, 
and is hardly practicable except for doublet, and 
perhaps for triplet levels. The difficulty lies in 
the fact that this procedure does not allow the 
states under consideration to be separated from 
all those of lower multiplicity. The simplest 
possible example is the triplet level of the He 
molecule. Here it is necessary to introduce two 
phantom orbitals, x and y, in addition to the real 
ones, a and b. There are then two structures, 


y x y—*, 
which must be considered. (It is not legitimate to 
neglect the second, as would seem reasonable at 
first sight.) The secular equation is found in the 
usual way to be 


Q-W-ta 
Q-W+a 
or W=Q+a, Q-a. 


Here, as before, Q is the coulomb integral, and a 
is the single exchange integral between a and b. 
All exchange integrals involving the phantom 
orbitals x and y vanish. Of the two roots, one 
corresponds to the singlet, and the other to the 
triplet level, but, before we can tell which is 
which, we have to calculate the singlet energy 
independently. This is found to be W=Q+a, and 
consequently, the desired triplet energy is given 
by the second root, W=Q-—a. In the general 


. Cf. L. Pauling, J. Chem. Phys. 1, 280 (1933); L. Paul- 
ing and G, W. Wheland, ibid. 1, 362 (1933). 


case, the calculation can become quite compli- 
cated. In order to find the 20 quintet levels of an 
8-electron problem, for example, we would have 
to solve a secular equation of the 132nd degree. 
Of the 132 roots thus obtained, 28 would 
correspond to the 14 singlet states (each oc- 
curring twice), 84 would correspond to the 28 
triplet states (each occurring three times), and 
the remaining 20 would correspond to the desired 
quintet states. 

A more straightforward application of the 
Slater method to such cases consists in actually 
setting up the proper number of independent 
wave functions, in the manner given by Kimball 
and Eyring,® and then evaluating the various 
matrix elements which occur in the secular 
equation. This latter is quite a laborious process, 
however, since no simple rules have yet been 
worked out for the states of higher multiplicity, 
analogous to those given by Kimball and Eyring® 
and by Pauling!’ for the singlet states.* 

In contrast to the above, however, the method 
of spin valence is very simply applicable to 
systems of high multiplicity. The justification of 
the following procedure is quite  straight- 
forward, and will be left to the reader.'® We 
introduce a single phantom orbital L, no matter 
what the multiplicity of the molecule may be, 
and to it we bond all the unshared electrons. 
(This apparent violation of the Pauli exclusion 
principle is legitimate here in consequence of the 
purely formal nature of the treatment.) Thus, for 
the triplet state of the H2 molecule, we have the 
single structure b 


L 
The corresponding function is 


J 
x 


where /, and /_ are constants, which are trans- 
formed under a rotation of axes in the same 
manner as a@ and 8, respectively. The function V 
is readily seen to be a linear combination of the 
three components of the triplet. There are only 

“s Cf. M. Born, Ergebnisse d. Exakten Naturwiss. 1931, 
387. 
* Note added in proof: Since the above was written, a 


method has been given by R. S. Bear and H. Eyring, J. 
Chem. Phys. 3, 98 (1935). 
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TABLE IV. 


(bc) 


A+B 


A(Q—W)—-(A+C)a 


—B 


B(Q—W)+(A-2B+C)a 


C+B 


C(Q-—W) —-(A+C)a 


two permutations to be considered, the identity 
I, and the single interchange (ab), both of which 
leave the structure unaltered. The secular equa- 
tion is accordingly 


or 


Thus the energy is simply W=Q-—a, in agree- 
ment with the result above. 

A somewhat more complicated problem is that 
of the triplet states of cyclobutadiene. There are 
three structures to be considered 


|Q—W-a| =0. 


As before, a complete independent set of func- 
tions is obtained by arranging the orbitals in a 
circle and drawing all structures in which no 
two bonds intersect. The usual rules apply to the 
resolution of structures with crossed bonds, in 
terms of the canonical set. The functions are 
set up in a manner analogous to the one above 
for He. Corresponding to A, for example, we have 


1 


4(6)3 


[a(1)8(2) —B(1)a(2) J 


The calculation is given in Table IV. 
The corresponding secular equation is 
Q-W-a 0 —a 
a Q-—W-2a a =0. 
—a 0 Q-W-a 


The solutions are W=Q, Q—2a, Q—2a, in 


complete agreement with the results obtained by 


Hiickel.® 
As in the case of the singlet states, it would be 


possible to carry through a discussion of such 
matters as the reduction of the secular equation 
by taking advantage of the symmetry of the 
molecule, and the approximate methods of 
solving the secular equations. This would lead to 
nothing essentially new, however, and conse- 
quently need not be discussed here in any further 
detail. 
ORBITAL DEGENERACY 

The problem of orbital degeneracy can also be 
treated by the method of spin valence. As in the 
case of the Slater method, however, the treat- 
ment becomes too complicated to be of very 
much use, except in actual numerical calculations 
for quite simple systems. We can represent a 
general wave function as 


& 


(2%(2N)!)4 


where the w,’s are functions of the spatial 
coordinates (analogous to the u of Eq. (3)), the 
d's are constants to be chosen so as to give the 
best possible function x, and the other symbols 
have the same significance as before. Some 
simplification is afforded by the fact that 
whenever wu; assigns two electrons, say 7 and s, 
to the same orbital, then a; will vanish unless ¢ 
contains the factor [a(r)B(s)—B(r)a(s) ]. If we 
apply perturbation theory by setting (7—W)x 
=0 and by then proceeding in the same manner 
as before, we obtain 
=0 (l=1, 2, pw). 
ik j 
This set of simultaneous homogeneous linear 
equations for the a’s takes the place of the set 
(9), and, like it, gives rise to a corresponding 
secular equation from which the energy can be 
obtained. It is doubtful if this method offers any 
advantages over the usual one, with which it is 
equivalent. It is given here only for the sake of 
completeness. 
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This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing‘ Editor 


The Structure of Cyanuric Triazide 


Mr. E. W. Hughes has recently described a structure for 
cyanuric triazide,' which in general outline is similar to 
one which I found some time ago and which was briefly 
described by Sir William Bragg in Nature of last July.? It 
differs, however, from mine in certain important details. 
The atomic distances in the cyanuric ring are given as 
1.33A and all equal, thus indicating an oscillating double 
bond as in benzene. I find two distances 1.38A and 1.31A, 
corresponding to fixed single and double bonds alternating 
in the ring. In Mr. Hughes’ structure the azide group 
departs considerably from linearity, being bent 15° at the 
middle nitrogen atom, while in my structure the group is 
straight or very nearly so. In the azide group, the nitrogen 
atom, which is directly linked to the cyanuric ring is un- 
fortunately not completely resolvable in the Fourier pro- 
jection on the (0001) plane because of a similar nitrogen 
atom of the second molecule of the unit cell being separated 
from it in the plane of projection by only 0.92A, so that 
there is rather less precision in locating its center than is 
the case with-the rest of the atoms. For this reason, it is 
not possible to be certain of the exact linearity of the azide 
group, but even so I estimate that any departure from it 
could not exceed 3° or 4°. 

It is not possible to follow Mr. Hughes’ work very closely, 
since he has not published the figures which he used for 
the intensities in his Fourier analysis nor the numerical 
values which he obtained from the summations, nor does he 
give any indication of the closeness of the points in the 
(0001) plane at which he carried out his summations. 

I may point out, however, that Mr. Hughes’ estimates 
of intensity were visual only. My intensity measurements 
were made on a Robinson photometer® and in my analysis 
I used the intensities from 99 planes of the (0001) zone, 
whereas Mr. Hughes only used 66 values. The intensities, 
which he has omitted are not in all cases the weakest. 
Many correspond to quite strong reflections, some of them 
having intensities which I found equal to one-half of the 
highest value I obtained for any plane in the zone. The 
Fourier series, which he used, must therefore be very far 
from convergent and it has been our experience in this 
laboratory that unless the series is very nearly convergent, 
that is, unless none but the very weakest reflections have 
been omitted, we do not obtain well-defined or correct 
Positions for the atoms, nor do contour lines which follow 
faithfully the results of the calculations, these having been 
carried out at a sufficient number of points, appear smooth 
or circular round the atom centers. 
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not later than the 15th of the month preceding that of the 
issue in which the letter is to appear. No proof will be sent 
to the authors. The usual publication charge ($3.00 per 
page) will not be made and no reprints will be furnished free. 


My intensity measurements have recently been put on 
an absolute scale by a few accurate determinations using 
the ionization method and molybdenum rays, to eliminate 
absorption. This has revealed an interesting feature of the 
structure, namely, that the electron density at the centers 
of the atoms is unusually high, rising to above 8 and 10 
electrons per square Angstrom unit for carbon and nitro- 
gen, respectively. Now that the work is quite complete, I 
hope that it will soon be published in detail. 

I. Knaccs 

Davy Faraday Laboratory, 

The Royal Institution, 
Albemarle, St., 
London, W. 1, England, 
February 8, 1935. 


1 Hughes, J. Chem. Phys. 3, 1 (1935). 
2 Bragg, Nature 134, 138 (1934). 
3 Robinson, Proc. Roy. Soc. A130, 120 (1930). 


Raman Spectra of Fused Salts 


In order to assist in the interpretation of results obtained 
from a series of thermodynamic studies of fused salts, we 
have begun a search for Raman spectra in a number of 
salts in the fused condition. 

We report some preliminary measurements of carefully 
purified salts made with a spectrograph of comparatively 
low dispersion. The mercury arc radiation was filtered with 
either an aqueous solution of sodium nitrite or a solution of 
iodine in carbon tetrachloride. No photochemical decom- 
position of the fused zinc chloride was observed after con- 
tinued illumination using either filter. In the pure state 
these salts did not etch the Pyrex Raman cells even after 
being confined at over 500°C for several weeks. Exposure 
times varied from four to twenty-four hours. 

The 100 cm™ anti-Stokes line from the 4358A exciting 
line was not observed since it would lie too close to the 
4339A mercury line for resolution. The 285 cm frequency 
value is the least exact, only one comparatively poor plate 
being available for this frequency. The Stokes and anti- 
Stokes lines are of comparable intensity as might be ex- 
pected for the temperature used. 


Tem- Exciting 
peratures line Displaced frequencies (cm~!) 


(400°, 430° 4358A 23049; 9845 
and 500°C) 23348 anti-Stokes) 
4046A 285 (broad); 101+8 (Stokes) 
23145; 1545 (anti-Stokes) 
4358A 155+8; 79+10 (Stokes) 
137+20 (anti-Stokes) 


Substance 
ZnCle 


ZnBr: (440°C) 
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Krishnamurti! investigated the Raman spectra of a 
number of inorganic halides in the solid state at room 
temperature. Those substances with heteropolar binding 
showed no Raman spectra whereas homopolar substances 
such as the mercuric halides showed strong Raman lines. 
Zinc chloride and cadmium iodide, substances of inter- 
mediate polarity, whose crystal structure is of the layer 
lattice type showed weak Raman lines. Our data would 
classify zinc bromide in the intermediate group with zinc 
chloride and cadmium iodide, the higher electrical conduc- 
tivity of the zinc bromide being offset by its lower melting 
point. 

Krishnamurti! found a frequency shift for zinc chloride 
of 234 cm™, which is in satisfactory agreement with the 
higher frequency line found by us. The frequency shifts of 
297 and 396 cm“ found from concentrated aqueous solu- 
tions of zinc chloride by Hibben? are probably due to com- 
plex ions or molecules.*: 4 Furthermore, Braune and Engel- 
brecht® have shown that the displaced frequencies of the 
mercuric halides are nearly the same for the fused salt as 
for the salt at room temperature. Grassman® and Thatte 
and Ganesan’ obtained similar results for a number of 
nitrates. The fact that the Raman frequencies persist some 
hundred degrees above the melting point indicates that the 
frequencies found here are not due to vibrations in the 
crystal lattice but are rather due to oscillations in the zinc 
halide grouping. 

The results for both zinc chloride and zinc bromide do 
not indicate a linear type of molecule. The determination of 


the angle between the valence bonds and strength of bonds 
awaits additional data. 

We failed to observe Raman lines from fused LiBr, PbClo, 
PbBr2 and CdCle, but the experimental conditions for 
obtaining these lines were not as satisfactory as with the 
zinc salts. 


Epwarp J. SALSTROM 
Louis Harris 

Research Laboratory of Physical Chemistry, 

Massachusetts Institute of Technology, 
February, 27, 1935. 

1 Krishnamurti, Ind. J. Phys. 5, 113 (1930). 

2 Hibben, Proc. Nat. Acad. Sci. 18, 532 (1932). 

3 Damaschun, Zeits. f. physik. Chemie B16, 92 (1932). 

4 Braune and Engelbrecht, Zeits. f. physik. Chemie B11, 416 (1931). 

5 Braune and Engelbrecht, Zeits. f. physik. Chemie B19, 303 (1932). 


6 Grassman, Zeits. f. Physik 77, 616 (1932). 
7 Thatte and Ganesan, Ind. J. Phys. 8, 341 (1934). 


The Absolute Rate of a Chemical Reaction 


Various writers!2»* have attempted to formulate the 
factors determining the absolute’ rate of a chemical reac- 
tion. Toa first approximation the task is not difficult. There 
are two methods of approach. One method is to calculate 
the number of favorable collisions from kinetic theory. The 
other is to treat the colliding molecule as a complex and 
apply the laws of statistical mechanics and chemical equi- 
librium for what they are worth. The two methods corre- 
spond at every point. Thus, the heat of activation in the 
latter treatment becomes the factor which determines the 
number of collisions of sufficient violence in the former 
method. The collision radius in the kinetic theory deter- 
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mines the moment of inertia and rotational entropy of the 
complex in the thermodynamic treatment. The greater the 
entropy of the complex, the greater must be the probability 
of reaction. By an arbitrary choice of collision radii, the 
“steric’’ factor may often be made unity, but this does not 
really give us any more information than we had before. 
Both methods involve uncertainties due to the lack of 
exact knowledge of force constants, collision radii, etc. 

Recently Eyring* has proposed to treat the subject by a 
method which is essentially the second method mentioned 
above. His method is not, therefore, novel, nor does it 
produce any results new to statistical mechanics. What he 
does propose that is new is to calculate the contours of the 
potential surface by wave mechanics so that the various 
uncertainties mentioned above may disappear. It may even 
be possible to calculate the rate of reflection at a barrier. 
This latter quantity may be described as a true steric 
factor. 

One may await results to see how well it will succeed. 
Presumably no one would claim that it is now possible to 
calculate the energy of any but the simplest chemical 
bonds with any precision by wave mechanics. Certainly a 
survey of the amount of work required to reach a satis- 
factory value for the heat of dissociation of the hydrogen 
molecule does not inspire one with confidence. The method 
of Heitler and London may be a better approximation in 
the case of a reactive complex than a stable molecule; it 
would be very fortunate if this is the case. It is, of course, 
possible to adapt calculations to experimental data, but in 
doing this one may produce agreement by changing one 
factor when some other factor is really at fault. 

Nevertheless, the writer has believed for a long time that 
it should be possible to determine all the ‘constants of a 
reaction rate from measurements at a single temperature. 

W. H. RopEBusH 

University of Illinois, 

February 27, 1935. 

1 W. H. Rodebush, J. Am. Chem. Soc. 45, 606 (1923); J. Chem. Phys. 
1, 440 (1933). 

439.(1928). and E. Wigner, Zeits. f. physik. Chemie (Haber band), 


30. K. Rice and H. Gershinowitz, J. Chem. Phys. 2, 853 (1934). 
4 Eyring, J. Chem. Phys. 3, 107 (1935). 


The Raman Spectrum of Deuterium 


A number of observers have studied the emission spectra 
of mixtures of hydrogen and deuterium. The Fulcher bands 
for HD and Dz have been photographed.! In the extreme 
ultraviolet Jeppesen,? and Beutler and Mie,* have obtained 
emission spectra from hydrogen-deuterium mixtures and 
have ascribed several bands to HD. From these data on 
HD, Urey and Teal‘ have calculated molecular constants 
for the molecule Dz in the ground state. No experimental 
data for this molecule in this state appear to have been 
obtained. 

We have obtained Raman spectra of gaseous deuterium 
(99.5 percent Dz) at a pressure of 2.5 atmos., using Hg 2537 
as an exciting line. 24- to 48-hour exposures were made. 
Two rotation lines and the Q branch of the vibration line 
were obtained. Their frequencies and mean deviations from 
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TABLE I. Observed and calculated Raman frequencies of D2. 


v obs. v calc. 


Difference 


179.6+0.5 
298.3 1.3 


2989.5 +1.3 


means in measuring plates are given in Table I, together 
with frequencies calculated from the equation of Urey and 
Teal, for the rotation-vibration energy levels of the mole- 
cule in the ground state. Using a formula given by Placzek® 
and due to Mannenbeck, we have calculated the relative 
intensities J of four of the strongest lines in the vibrational 
Q branch. The result gives 2990.0 cm~ for the position of 
the center of gravity of the Q branch, and this is in agree- 
ment with the observed frequency. 

The constants of Urey and Teal are therefore confirmed. 
This result is important in connection with the use of 
deuterium to establish the structure of hydrogen-contain- 
ing molecules.* In such cases M—D distances are assumed 
equal to M—H distances, where M is any atom. The 
greatest deviation from this assumption would probably 
occur in Do; and since no deviation is found in this case the 
assumption appears to be justified. Moreover the force 
constant for D2 is the same as that for He, and this suggests 
that M—D force constants are equal to those for M—H. 

Tuomas F. ANDERSON 
Don M. Yost 
California Institute of Technology, 
March 2, 1935. 

1M. F. Ashley, Phys. Rev. 43, 770 (1930); G. H. Dieke and R. W. 
Blue, oe 133, 611 (1934); Phys. Rev. 47, a ¢ 935). 

R. Jeppesen, Phys. Rev. 45, 480 (1934); 44, 165 (1933). 

2H. Beutler and K. — Naturwiss. 22, 418 (1934); K. Mie, Zeits. f. 
ony 91, 475 (1934); H. Beutler, Zeits. f. physik. Chemie B27, 287 

‘H. C. Urey and G. K. Teal, Rev. Mod. Phys. 7, 34, (1935). Table 
XV, p. 80. A complete set of references is given in this review. 

5G. Placzek, Rayleigh-Streuung und Raman-Effekt (Leipzig, 1934) 
Eqs. (18), (15a), p. 343. 

F. Rasetti reports a line v10% for hydrogen. Phys. Rev. 34, 367, 
(1929). The matter is worth further investigation. 


6 For example the structure of NHs: Barnes, Benedict and Lewis, 
Phys. Rev. 45, 347 (1934). 


A Transition in Potassium Superoxide 


During the study of the magnetic properties! of potas- 
sium superoxide, KO:, I noticed that on cooling in liquid 
air the substance changes from a deep orange to a light 
cream color. In order to see whether or not this change in 
color accompanies a transition from one form to another, 
and if so to determine the transition temperature, the 
following experiments were carried out. 

Samples of the superoxide and of sodium chloride were 
put in a glass tube 7 mm in diameter, each in a layer 2 cm 
deep, with the two junctions of a differential thermocouple 
at the centers of the two layers, and the junction of another 
thermocouple (for measuring the temperature) at the center 
of one of the layers. The tube was sealed from the atmos- 
Phere by a plug of cotton followed by a layer of picien. 
The tube was enclosed in a copper jacket about 5 mm thick, 
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in order to insure a uniform heat transfer between the cell 
constants and the surroundings. Mercury was added as a 
further aid to this end and the whole was enclosed in a test 
tube. The thermocouples were single junction and were 
made of No. 32 copper and constantan wire joined with 
silver solder. The differential thermocouple was connected 
directly to a Leeds and Northrup Type R galvanometer 
which had a period of 5.5 sec. and a sensitivity of 2.7 mm/m 
per mv. The galvanometer scale was 60 cm from the instru- 
ment. No damping resistance was found necessary with this 
instrument. This differential thermocouple was used only 
as a means of observing differences of temperature between 
the oxide and the chloride. The actual temperature of the 
oxide was measured by means of the other thermocouple 
which was used in conjunction with a Leeds and Northrup 
Type K potentiometer. This thermocouple was calibrated 
by comparison with nine fixed points. 

Warming curves were obtained from liquid air tempera- 
ture to the temperature of melting ice. The cell was first 
cooled in liquid air for about an hour, the air was then 
removed from the enclosing Dewar flask, a cotton plug 
was inserted into the mouth of the flask to prevent irregular 
heating due to convection, and the cell was allowed to rise 
to room temperature. The most uniform curves were ob- 
tained when the time of warming over the temperature 
range studied was between 4 and 5 hours. 

Cooling curves were obtained by enclosing the cell in an 
unsilvered Dewar flask and cooling by immersing this flask 
in liquid air. The flask that was used allowed the tempera- 
ture of the cell to fall from room temperature to liquid air 
temperature in about 4 hours. 

The warming curves and cooling curves all showed a pro- 
nounced break, corresponding to a transition at the tem- 
perature —75.5°+0.5°C. (No break in this region was 
shown by control curves obtained by replacing the super- 
oxide in the cell by sodium chloride.) Observation of the 
color change showed that it occurs at the same tempera- 
ture, to within about +15°. 

The experimental method is not sufficiently delicate to 
distinguish between a true polymorphic transition and a 
gradual transition;? if the transition is of the latter type, 
it is reasonably sharp, extending over a temperature range 
not greater than three or four degrees. 

The normal superoxide ion is in the state °II, and the 
separation of the two levels *IIj/2 and °Iz/2 is probably 
about the same as for NO, 121 cm~. It is interesting to 
note that this energy value is equal to kT at the tempera- 
ture T=213°A, only 15° from the observed transition tem- 
perature; though until more definite evidence is obtained 
the suggestion that the transition of the solid substance is 
related to the electronic transition of the molecules must be 
considered as speculative. 

The author is grateful to Professor Linus Pauling for his 
interest in the investigation. 

E. W. NEuMAN?® 


California Institute of Technology, 
March 15, 1935. 


1E. W. Neuman, J. Chem. 
2 Linus Pauling, Phys. Rev. 36, 4 


3 Now at Carroll College, Capitol iL Helena, Montana. 
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